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CHAPTER  1 
INTRODUCTION 

Energy  conservation  is  an  important  item  for  consideration  to  all, 
building  designers,  developers,  owners  and  managers.     Air-to-air  heat  re- 
covery, through  processed  air  with  energy  recovered  from  exhaust  air,  is 
one  area  of  investigation  for  potential  energy  and  cost  savings. 

Although  ideally  suited  for  hospitals,  office  buildings,  apartments, 
hotels,  schools,  auditoriums  and  factories,  air-to-air  heat  recovery  equip- 
ment is  applicable  to  almost  any  commercial,   institutional,  or  industrial 
facility.     Therefore,  it  should  be  carefully  considered  wherever  wasted 
energy  is  discharged  from  a  process  or  building. 

While  the  reserve  of  energy  sources  is  not  unlimited,  it  is  being  wasted 
in  enormous  quantities. 

But  wherever  energy  is  wasted,  there  is  an  opportunity  to  recover 
some  portion  for  some  useful   purposes.     This  was  the  theme  of  a  paper 
presented  by  P.  C.  Greiner,  Edison  Electric  Institute,  New  York  at  the 
recent  National   Electric  Heating  Conference  held  in  Cincinnati   [1]*.     Very 
simply,  he  says,  "Don't  buy  energy  twice". 

Heat-loss  through  air  change  represents,  in  general,  some  30  percent 
of  the  overall   heat  loss  in  a  multiple-dwelling  building  and  approximately 
20  percent  of  such  loss  in  a  single  family  house  [2]. 


Number  in  brackets  designate  references  in  the  list  of  References. 


VARIOUS  HEAT-RECOVERY  SYSTEMS  CURRENTLY  AVAILABLE  FOR  DOMESTIC  USE: 

Existing  heat  recovery  systems  are  surveyed  by  examining  the  advan- 
tages and  disadvantages  of  each  from  various  points  of  view,  particularly: 

-  that  of  efficiency,  which  is  obviously  of  capital   importance; 

-  that  of  pressure  loss,  which  increases  ventilation  costs  and  make  acoustic 
problems  more  difficult  to  resolve; 

-  that  of  frost  occurrence,  which  affects  efficiency; 

-  that  of  sanitation:   the  forced-air  intake  must  not  be  polluted  by  the 
exhausted  air. 

1.     Systems  Based  on  the  Thermal   Capacity  of  a  Material: 

Figure  (la)  represents  a  rotary  system,  the  most  typical   type  of 
these  systems.     Recirculation  of  a  volume  of  gas  occurs  at  the  time  of  air- 
current  inversion  or  reversal.     This  disadvantage  can  be  minimized  by  purging 
the  heat  recovery  system  at  the  time  of  reversal,  employing  some  additional 
apparatus.     In  addition,  the  heat  recovery  system  may  gradually  become 
fouled  and  clogged.     This  problem  can  be  overcome  either  by  using  a  washing 
device  or  by  fitting  filters;  however,  both  solutions  have  some  disadvantages. 

The  rotary  heat  recovery  system  consists  of  a  porous  wheel  which 
rotates  about  an  axis  parallel   to  two  gaseous  streams.     The  streams  are 
separated  by  a  wall   that  diametrically  intersects  the  wheel.     With  this 
apparatus,  the  air  streams  follow  a  fixed  course,  and  the  collector 
(accumulating)  material  passes  from  one  stream  to  another  at  variable  speed. 
This  system  under  proper  conditions  is  80  percent  efficient  [2]  and  has  the 
same  disadvantages  as  the  stationary  porous  system.     A  purger  permits 
limiting  the  recirculation  of  stale  air  at  the  point  where  the  wheel   crosses 
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Figure  1.     Schemes  of  Several   Principles  of  Air-to-Air  neat  Exchangers. 


the  boundary  within  the  stream  of  exhausted  air  and  fresh  air.     But  there 
is  a  risk  of  recirculation  on  a  level  with  the  seal  between  the  rotating 
part  and  the  barrier  separating  the  two  streams.     This  risk  can  usually  be 
limited  if  the  pressure  of  the  fresh  air  exceeds  the  exhausted  air. 

If  the  moist  exhaust  air  is  below  32°F,  as  it  passes  over  that  part  of 
the  system  material  whose  temperature  is  close  to  the  outside  temperature, 
frost  is  apt  to  form  within  these  heat  recovery  systems.     This  condition 
can  be  improved  by  use  of  a  pre-heating  system  or  by  shutting  the  system 
off  below  a  certain  outside  temperature.     Unfortunately,  these  remedial 
steps  decrease  the  systems  efficiency. 

Typical  recovery  efficiencies  for  non-hygroscopic  rotary  are  70-80 
percent  sensible  heat  [3]. 

2.     Systems  using  Direct  Heat  Transfer: 

Figure  (lb).     These  type  of  exchangers  may  be  of  metal,  plastic,  or 
combined  fiber  materials.     Air  flow  is  usually  either  cross  or  counter 
flow.     Counter  flow  provides  the  greatest  temperature  difference  for  maximum 
heat  transfer,  and  cross  flow  can  allow  more  convenient  air  connections. 
The  usual  plate  exchangers  transfer  only  sensible  heat,  when  the  temperature 
of  one  air  stream  is  lower  than  the  dew  point  of  the  other,  and  direct  conden- 
sation occurs.     Preheating  to  avoid  condensation  could  prevent  unusually 
high  pressure  drop,  corrosion,  or  freeze-up  due  to  presence  of  water  on  the 
surfaces.     As  desired,  recovery  may  be  limited  and  controlled  by  having  a 
portion  of  the  supply  air  by-pass  the  exchanger  by  means  of  dampers.     There 
is  no  leakage  between  exhaust  and  supply  air  if  air  passages  are  sealed. 
Usual  recovery  efficiencies  are  40  to  60  percent  sensible  heat  [3]. 


3.  Systems  Employing  an  Intermediary  Energy-Carrying  Fluid: 

(Run-Around  System) 

Figure  (lc).     This  system  consists  of  an  energy-carrying  fluid  circu- 
lating between  two  finned  tube  coils,  one  set  within  the  warm-air  stream, 
the  other  within  the  cold  air  stream.     A  pump  serves  to  circulate  this  fluid. 

This  system  offers  significant  advantage  over  the  previous  one: 

-  exhaust  ducts  and  forced-air  intake  ducts  may  be  placed  at  some  distance 
from  each  other, 

-  contamination  of  bad  odors  is  not  likely  to  occur. 

Frost  problems  and  fouling  difficulties  are  resolved  as  in  the  first 
system. 

4.  Run-Around,  Multiple  Tower,   Desiccant  Spray  Exchanger: 

Figure  (Id).     This  system  is  similar  to  the  run-around,  coil  and  pump 
system.     An  extended  packed  surface  is  substituted  for  the  cooling  coils, 
such  as  cooling  tower  fill.     A  liquid  absorbent  (usually  a  solution  of  the 
halogen  group  such  as  lithium  chloride/water)  is  substituted  for  the  water/ 
anti-freeze  circulating  liquid.     Another  difference  is  that  the  absorbent 
liquid  is  sprayed  counter  flow  to  the  air  streams  through  the  extended 
packed  surface.     This  necessitates  two  solution  pumps  to  complete  the  run- 
around  circuit. 

This  system  provides  total  heat  transfer  as  the  solution  absorbs   (or 
desorbs)  heat  and  water  vapor  from  the  exhaust  air  stream.     Moreover,  by 
circulation  and  contact  the  circulating  liquid  with  the  system,  the  supply 
air  stream  is  preconditioned  to  approach  the  system's  temperature  and 
vapor-pressure  difference.     The  system  is  also  reversible  because  of  its 


nature  of  action,  similar  to  the  hygroscopic  rotary  recovery  system.  It 
will  pre-cool  and  dehumidify  in  the  summer  and  will  preheat  and  humidify 
in  the  winter. 

5.     Heat  Pipe  Air-to-Air  Heat  Transfer  System: 

Figure  (le).     In  this  system,  the  ducts  are  arranged  such  that  inlet 
air  stream  flows  through  one  side  of  the  unit,  and  the  outlet  air  stream  flows 
through  the  opposite  side.     Heat  is  thus  transferred  from  one  stream  to 
another  stream. 

The  unit  is  constructed  of  a  series  finned  tubes.     Each  tube  consists 
of  an  envelop  (the  tube),  a  wick,  and  a  working  fluid  and  is  sealed  at  both 
ends.     These  tubes  are  the  actual   heat  pipes.     Heat  applied  to  one  end  of 
the  unit  evaporates  the  working  fluid  from  the  wick.     The  vapor  then  flows 
to  the  cold  end  of  the  tube,  where  it  is  condensed  and  returned  by  the 
wick  to  the  hot  end  for  re-evaporation,  thus  completing  the  cycle. 

The  heat  pipe  is  a  completely  reversible,  essentially  iso-thermal 
device.     When  a  section  of  the  envelop  is  exposed  to  a  temperature  gradient 
with  respect  to  another  section,     the  cycle  commences  seizing  the  existing 
gradient. 
Advantages: 

1.  Manufacturers  claim  up  to  68%  efficiency  [4]. 

2.  No  cross  contamination  of  air  streams  occurs. 

3.  The  leaving  air  temperature  may  be  controlled  by  tilting  the 
coil . 


Disadvantages: 

1.  The  units  are  costly  [4]. 

2.  Side-by-side  air  discharge  and  intake  is  required,  which  necessi- 
tates expensive  duct  work  to  bring  the  supply  and  exhaust  air  to 
the  equipment  room  for  processing.     Some  disadvantage  exists  for 
all  but  run-around. 

3.  The  units  are  reportedly  limited  to  a  maximum  length  of  8  ft  [4]. 

4.  A  preheating  ceil    is  required  for  temperatures  below  32     F  outside 
temperature  to  prevent  the  exhaust  air  condensate  from  freezing. 

In  this  study,  the  purpose  is  to  develop  a  technical  and  economic 
analysis  of  the  Run-Around  Coil  System  and  Heat  Pipe  System  heat  exchanger 
types  used  in  exhaust  air  heat  recovery. 

These  two  are  typical   of  all   types  available.     The  analysis  of  the 
heat  pipe  system  depends  on  performance  data  reported  by  manufacturer 
except  for  the  "run-around"  system.     The  "run  around"  system  is  a  "built-up" 
system  of  standard  components  available  from  different  manufacturers. 


CHAPTER  2 

"RUN-AROUND  COIL  SYSTEM" 
General   Description 

Figure  (lc).     With  the  standard  extended  surface,  finned-tube  water 
coils  may  be  used  with  a  circulation  pump  for  recovery  of  sensible  heat 
from  exhaust  air.     One  coil   is  located  in  the  exhaust  air  stream  and  the 
other  in  the  supply  air  stream.     A  pump  circulates  water  or  a  water-antifreeze 
solution  between  the  coils,  transferring  energy  from  the  exhaust  recovery 
coil   to  the  supply  air  coil.     The  direction  of  the  flow  of  solution  to  the 
air  is  usually  counter  for  maximum  heat  transfer. 

The  system  is  seasonally  reversible  (i.e.,  the  air  to  be  pre-conditioned 
by  the  supply  air  coil  will   reflect  an  approach  to  the  solution  temperature 
from  the  exhaust  air  coil,  thereby  preheating  when  the  outdoor  air  is  cooler 
than  the  exhaust  and  precooling  when  the  outdoor  air  is  warmer  than  the 
exhaust) . 

This  system  operates  generally  for  sensible  heat  recovery  only,  except 
when  the  temperature  of  one  air  stream  is  lower  than  the  dew  point  temper- 
ature of  the  other,  and  direct  condensation  occurs.     Preheating  the  lower 
air  stream  to  avoid  condensation  will  prevent  unusually  high  pressure  drop, 
corrosion,  or  freeze-up  due  to  water  on  the  surfaces.     When  required, 
recovery  may  be  limited  and  controlled  through  the  use  of  face  and 
by-pass  dampers.     Components  are  usually  selected  for  ratings  that  will 
achieve  recovery  efficiencies  of  40  to  60  %  sensible  heat  [3]. 


Typical   Applications 


Comfort 

Offices 

Laboratories 

Schools 

Plants 

Recreational   Facilities 

hospitals 
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Restaurants 


Industrial /Process 
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Theory  ar.d  Modeling  Equations 

Heat_ aajarca  and  Transfer  Sanations 

The  heat  transfer  ana  balance  equations  for  this  system  were  mrltter 
as  (see  notation  of  Figure  2): 


Coil   1 


Coil   2 


Exhaust  air 


Figure  (2) 
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(1) 


(T5-T2)  -    (T?-T4) 
«2  •   Cpa   (W  =  U2A2  «i[(T5-T2)/(T7-l4)J 

(t5-t3)  -  (yy 

=  U1A1    £n[(T5-T3)/(T6-l1)J 

=  Q.    •  C      (To-T,)  (Winter  operation)         (2) 

assuming  counter- flow  operation 

Q3  +  q5  -  Q4  (3) 

G3T6  +  G5T5  =  G4T7  (4) 

By  trial-and-error,  given  values  of  T1 ,  Tg,  G1 ,  G2>  U1  and  U2>  a 
balance  case  can  be  obtained  and  Tg,  Tg,  T?,  T3  and  T4  can  be  determined. 
Details  for  their  determination  are  given  in  Appendix  A. 

Frost  Control   Determination 

Frost  control   is  needed  only  during  winter  operation.     If  T?  <  32    F, 
Q5  may  be  selected  to  require  T?  ■  32  °F.     More  details  are  given  in 
Appendix  A. 

Make  up  and  Exhaust  Air  Coils 

The    U     value  of  the  coil  based  on  outside  (air-side)  area: 

u       VA°-58  166  VEG0-8 


The  air  pressure  drop: 


,_2  ,  n    0.4    ,    0.6 

^-\ fa  ft  (6) 


APD  =  27.69  x  *— fg f-  c    , 

(5.22xlOlu)   DvS*s   ^V         ^V 

The  pressure  drop  of  the  ethylene  glycol: 

PEG  =  1.495  (^p)  [0.0039  WZ  +  0.0875  (W-l )  +  0.3]  + 

35.78^  (7) 

Equation   (5)  is  that  developed  by  Stoecker  [5].     Details  of  determin- 
ation of  equations  (6)  and  (7)  are  given  in  Appendix  A. 

Ethylene  Glycol 

The  mass  flow  rate  and  heat  capacity  are: 

Q  ■  176431.84  VEG  DI2  (1+0. IX)  (2) 

G  =  (1-0.398  X)  Q  (9) 

Details  of  development  of  equations  (8)  and  (9)  are  given  in 
Appendix  A. 

Energy  Recovered 

The  system  is  built  for  sensible  heat  transfer  in  which  no  moisture 
is  exchanged  between  air  streams.  The  sensible  heat  gain  by  the  outside 
air  is: 

q  ■  1.08  cfms   Oyty  (10) 


Regarding  the  sensible-heat  recovered  per  season,  information  from 
figures  (3),  (4),  (5)  and  (6),  which  gives  temperature  frequency  during 
heating  and  cooling  seasons  for  an  area  averaging  4000  to  5000  and  5500 
to  6500  degree  days  respectively,  have  been  used. 
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Operating  Energy 

The  energy  required  to  operate  the  fans  and  pump  is 


Fan: 


0.4003876  cfmxAPD 

1  = iff. 


(11) 


Pump: 


0.642676  GPU  x  S  x  PEG 

i = m: 


Development  of  equations  (11)  and  (12)  are  in  Appendix  A. 


(12) 


Net  Energy  Recovered 

q  =  q  (energy  recovered)  -  q  (operating  energy  for 
fans  and  pump) 


(13) 


CHAPTER  3 

THE  COST  EQUATIONS 

To  determine: 

The  Initial   Cost  for 

Coil 

4.   it,       **   \  19.1   -   0.083  Z  ,-,.) 

cost     S/sq.ft.     =  H71 0~73  l     ' 

Wl/3  (A0/1000)U-U 

Estimated  cost  for  standard  air-coil   exchanger  in  $/sq.ft.  of  bare-tube 

surface  as  formulated  by  Edward   [8]. 

The  Total    Initial   Cost 

cost  ($)  =  I  cost  (coils  +  pump  +  labor  +  ethylene  glycol  +  fans) 

(15) 
The  Operating  Energy  Cost  for 
Fan 

cost  ($/year)  =  (Fan  BHP)  x  0.746  gj|r  x  -^3—  (Elec.  cost)  x  yjjjjf 

(16) 
Pump 

cost  ($/year)  «  (Pump  BHP)  x  ""^p^  x  ^$   (Elec.  Cost)  x  ^ 

(17) 
The  Total  Operating  Energy  Cost 

cost  ($/year)  =  Fans  operating  energy  cost  ($/year)  + 

Pump  operating  energy  cost  ($/year)      (18) 
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The  Capital  Equipment  Reduction  for 
Boiler 


rost  m  =  Btuh  energy  recovered  in  winter   ^_  Average  equipment  cost 
cost  i»J    100000  Btuh/Therm  Therm 

(19) 
Chiller 

cost  ($)  =  Btuh  energy  recovered  in  summer   J_  A     equipment  cost 
cost  i$;     120Q0  Btun/Ton  Ton     M   M 

(20) 

The  Total  Capital  Equipment  Reduction: 

cost  ($)  '   $  Boiler  +  $  Chiller  (21) 

The  Energy  Gross  Recovery: 

Winter 

,&l         ,   Btuh  energy  recovered  in  winter    $   Hn.itinn  rnst 
savings  ($/year)  =   100ooo  Btuh/Therm x  Therm  Heatlng  Cost 

(22) 

Summer 

it i         \       Btuh  energy  recovered  in  summer      _$_  r„-n„n  rn.t 
savings  ($/year)  =      12000  B^n/Ton *  To7  Coolln9  Cost 

(23) 
The  Total   Energy  Gross  Recovery: 
savings  ($/year)  =  Winter  G.R.   ($/year)  +  Summer  G.R.   ($/year)         (24) 

The  Net  Energy  Gross  Recovery: 

savings  ($/year)  ■  Total   G.R.   ($/year)  -  Total  operating  ($/year) 

(25) 

The  Years  of  Pay  Back: 

_     .  ,    ,  %   Total  Initial  Cost  -  $  Total  Capital  Reduction 
Pay  back  (yrs.)  -  Net  Gross  Recovery  ($/Yrs)       ^ 

Details  of  equations  (14  to  26)  are  given  in  Appendix  A. 
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Calculational  Scheme 

A  computer  program  (shown  in  Appendix  B)  was  written  which,  by  in- 
putting the  coil   descriptions,  air  flow  rates,  air  temperatures,  fan 
efficiencies,  pump  efficiency,   initial  system's  component  cost,  unit 
heating  and  cooling  energy  cost,  and  an  optional  water  antifreeze  solution 
percentage  and  velocity,  allows  the  design  engineer  to  estimate  the  net 
energy  recovered  and  the  pay  back  period.     The  program  iterates  between 
the  make  up  and  exhaust  coils  of  the  system  until   the  antifreeze  solution 
temperatures  (Tg,  Tg  and  Ty)  balance.     At  this  condition  all   energy  and 
dollar  quantities  are  calculated. 

The  program  determines  the  specifications  for  controlling  condensate 
frosting  at  the  design  outside  air  temperature  and  calculates  the  resulting 
capacity  with  frost  control  operating. 


CHAPTER  4 

SPECIFIC  SYSTEM  ANALYSIS 
Results  and  Discussion 

The  potential  savings  can  be  evaluated  for  a  heat  recovery  system. 
An  optimum  system  consists  of  the  most  favorable  combination  of  the  following 
variables: 

1 .  Length  of  the  coil ,  Z. 

2.  Height  (or  number  of  tubes  high)  of  coil,  HT. 

3.  Number  of  rows  of  tubes  deep  (parallel  to  the  path  of  air  flow),  W. 

4.  Fin  spacing,  FN. 

5.  Glycol   flow  rate,  VEG. 

The  qualitative  effects  of  the  previously  mentioned  five  variables 
are  as  fol lows : 

Increasing  the  length    or   number  of  tubes  high,  increases  the  heat 
transfer  area  but  also  increases  the  cost.     Furthermore,  increasing  the 
length    or   the  number  of  tubes  high  not  only  increases  the  cross-sectional 
area  for  air  flow,  reducing  the  velocity  and  decreasing  the  air-side  heat 
transfer  coefficient,  but  also  decreases  the  power  required  for  the  fan. 
The  first  and  second  above  variables,  were  combined  as  a  function  of  one 
variable,  air  face  velocity  (VA)  as  shown: 

VA  =  -£EM 

VH      HTxZ 

Increasing  the  number  of  rows  of  tubes  deep  increases  the  heat  transfer 

area,  but  also  increases  the  initial  cost  of  the  coil  as  well  as  the  cost 

of  pumping  both  air  and  ethylene  glycol. 
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Spacing  the  fins  closer  together  increases  the  cost  of  the  coil  and 
the  air  pressure  drop,  but  also  increases  the  heat-transfer  area. 

Finally,  a  high  flow  rate  of  the  ethylene  glycol  increases  the  glycol 
side  heat-transfer  coefficient,  but  also  increases  the  pumping  cost. 

Values  of  fin  spacing,  face  velocity,  ethylene  glycol  velocity  and  coil 
depth,  were  selected  for  the  application  considered  here  to  be  10  fins  per 
inch,  500  fpm,  8  fps  and  4  rows  respectively. 

These  four  parameters  were  chosen  to  study  the  individual  effects  on 
the  system  by  subjecting  each  parameter  to  a  variation  while  keeping  the 
other  parameters  at  their  standard  values. 

The  effect  of  the  fin  spacing  (FN)  variation  on  annual  energy  recovered, 
net  annual  energy  recovered  cost  and  pay  back  periods,  are  shown  in  figures 
7,  8  and  9  respectively.     The  optimum  and/or  economic  value  of  fin  spacing 
in  all  these  figures  is  10  fins  per  inch. 

In  figure  7,  total   recovery  unit  area  (A)  curve  shows  a  gradual   in- 
crease as  the  number  of  fins  per  inch  increases.     The  reason  is  that  the 
total   recovery  unit  area  is  directly  proportional  to  the  number  of  fins  per 
inch. 

Total  annual  operating  energy  (TAOE)  curve  shows  a  gradual   increase 
as  the  number  of  fins  per  inch  increases.     The  reason  is  that  total  annual 
operating  energy  is  directly  proportional   to  air  pressure  drop  (APD)  and  the 
air  pressure  drop  is  directly  proportional   to  the  number  of  fins  per  inch. 

Total  annual   gross  recovery  (TAGR)  curve  shows  a  gradual  increase  until 
it  reaches  the  optimum  value  of  the  number  of  fins  per  inch.     Further  in- 
creases in  the  number  of  fins  per  inch  shows  a  small   response.     The  reason 
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is  that,  total  annual  gross  recovery  is  a  function  of  overall   heat  trans- 
fer coefficient  (U)  and  total   recovery  unit  area.     Overall  heat  transfer 
coefficient  is  inversely  proportional   to  the  number  of  fins  per  inch  as 
shown  in  equation  (5).     The  rate  of  decreases  in  overall  heat  transfer 
coefficient  is  proportionally  much  smaller  than  the  increase  in  total  re- 
covery unit  area.     However,  the  numerical   results  show  that  the  by-pass 
was  open  for  preventation  of  frost  formation.     This  limited  total  annual 
gross  recovery,  as  noted  in  the  figure  by  reducing  the  ethylene  flow  through 
the  make  up  coil . 

Total  annual  net  gross  recovery  (TANGR)  is  the  difference  between  the 
total   annual  gross  recovery  and  the  total  operating  energy,  as  is  shown  in 
the  curve. 

Figure  8,  shows  the  effect  of  annual  net  energy  recovered  cost  for  the 
oil,  gas,  and  electric  system  (if  it  is  used  for  heating)  with  respect  to 
the  variation  of  the  fin  spacing. 

The  annual  net  energy  cost  recovered  is  directly  proportional   to  the 
total   energy  gross  recovered  cost  (TGR)  and  inversely  proportional   to  total 
operating  energy  cost  (TOE).     Total   energy  gross  recovered  cost  for  oil, 
gas,  and  electric  are  functions  of  the  total  annual  gross  recovery  and 
heating  costs  of  oil,  gas,  and  electric.     The  total  operating  energy  cost 
is  a  function  of  the  total  annual  operating  energy  and  the  electric  operating 
costs . 

The  annual  net  energy  recovered  cost  for  oil  curve  shows  an  increase 
until  the  number  of  fins  per  inch  value  equal  to  10.  Further  increase  in 
the  number  of  fins  per  inch  value  results  in  an  appreciable  drop  in  the 
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annual  net  energy  recovered  cost  for  oil.     The  behavior  of  the  curve  is  the 
difference  between  the  total  annual   gross  recovered  multiplied  by  oil  heating 
cost  (taken  0.26  $/Therm)  and  the  total  annual  operating  energy  multiplied 
by  electric  operating  cost  (taken  0.025  $/KW-HR). 

There  is  an  overall  decrease  in  the  annual  net  energy  cost  recovered 
curve  by  increasing  the  number  of  fins  per  inch.     However,  up  to  the  fin 
spacing  value  equal   to  10,  the  decrease  is  more  gradual   than  for  higher 
values  of  the  fin  spacing.     The  behavior  shown  in  the  curve  is  due  to  low 
gas  heating  cost  (taken  0.13  $/Therm). 

There  is  a  noticeable  increase  in  the  value  of  the  annual  net  energy 
cost  recovered  for  electric  up  to  the  fin  spacing  equal  to  10.     For  values 
of  the  fin  spacing  greater  than  10,  a  very  small  increase  is  seen  as  evident 
from  the  curve.     That  is  due  to  hign  electric  heating  cost  (taken  0.73 
$/Therm). 

The  noticeable  difference  in  the  slopes  of  the  annual  net  energy  re- 
covered cost  for  oil,  gas  and  electric  curves,  with  respect  to  the  vari- 
ation of  the  fin  spacing,  are  due  to  lowest  heating  cost  for  gas  and  oil, 
and  the  highest  heating  cost  for  electricity. 

The  pay  back  period  is  directly  dependent  on  initial  costs   (EIC),  and 
is  inversely  proportional   to  total   capital  equipment  reduction  (TECR)  and 
net  energy  recovered  costs.     The  pay  back  for  oil,  gas,  and  electric  (PAYB0, 
PAYBG,  and  PAYBE)  curves  are  represented  in  figure  9.     As  is  evident  from 
these  curves,  the  minimum  pay  back  period  occurs  at  the  fin  spacing  value 
of  10  fins  per  inch. 
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Figures  10  and  11  represent  the  relation  between  the  total  annual  gross 
recovered,  total  annual  net  gross  recovered,  total  recovery  unit  area,  total 
annual  operating  energy,  the  annual  net  energy  cost  recovered  for  oil,  gas 
and  electric  against  the  face  velocity  (VA). 

The  total  recoveryunit  area  decreases  as  the  face  velocity  increases, 
as  is  evident  from  the  graph.  The  total  recovery  unit  area  is  directly  pro- 
portional to  the  length  of  the  coil.  The  face  velocity  is  inversely  pro- 
portional to  the  product  of  the  height  and  the  length  of  the  coil.  Keeping 
the  height  constant,  the  length  decreases  as  the  face  velocity  increases, 
hence  the  above  result  between  the  total  recovery  unit  area  and  the  face 
velocity. 

The  total  annual  operating  energy  increases  as  the  face  velocity 
increases.  The  total  annual  operating  energy  is  directly  proportional  to  the 
pressure  drop.  The  face  velocity  is  directly  proportional  to  the  pressure 
drop.  This  makes  a  noticeable  increase  in  the  total  annual  operating  energy 
as  the  face  velocity  increases. 

The  total  annual  gross  recovery  decreases  as  the  face  velocity  increases. 
The  total  annual  gross  recovery  is  directly  proportional  to  the  product  of 
the  overall  heat  transfer  coefficient  and  the  total  recovery  unit  area. 
The  overall  heat  transfer  coefficient  is  directly  proportional  to  the 
fractional  power  of  face  velocity  (VA   )  and  the  total  recovery  unit  area 
is  inversely  proportional  to  the  face  velocity.  Hence  UA  decreases  as  VA 
increases  and  thus  the  behavior  is  seen  in  the  curve  for  total  annual 
gross  recovery. 

The  total  annual  net  gross  recovered  decreases  as  the  face  velocity  in- 
creases. The  total  annual  net  gross  recovered  is  the  difference  between 
the  total  annual  gross  recovered  and  the  total  annual  operating  energy. 
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The  annual  net  energy  cost  recovered  for  oil,  gas  and  electric  curves 
decreases  as  the  face  velocity  increases.  The  annual  net  energy  cost 
recovered  is  directly  proportional  to  the  total  energy  gross  recovered  cost 
and  inversely  proportional  to  the  total  operating  energy  cost.  The  total 
energy  gross  recovered  cost  is  function  of  the  product  of  the  total  annual 
gross  recovered  and  heating  cost.  The  total  operating  energy  cost  is  the 
function  of  tne  product  of  the  total  annual  operating  energy  and  the  electric 
operating  costs. 

The  noticeable  difference  in  slopes  of  the  annual  net  energy  cost  re- 
covered for  oil,  gas  and  electric  is  due  to  different  oil,  gas,  and  electric 
heating  costs. 

From  the  above  curves,  it  is  seen  that  the  best  value  of  the  face  ve- 
locity is  between  400  to  600  fpm.  500  fpm  was  selected  for  a  standard  run. 

Figure  12  shows  that  the  total  recovery  unit  area  is  constant  as  the 
ethylene  glycol  velocity  increases  because  the  total  recovery  unit  area  is 
not  dependent  on  the  ethylene  glycol  velocity. 

The  total  annual  operating  energy  increase  as  the  ethylene  glycol  ve- 
locity increases.  The  reason  is  that  the  total  annual  operating  energy  is 
directly  proportional  to  the  pressure  drop  in  ethylene  glycol  (PEG).  The 
pressure  drop  in  ethylene  glycol  is  dependent  on  the  ethylene  glycol  velocity 
(VEG2  and  VEG1"75)  as  shown  in  equation  (7),  hence  the  above  result  between 
the  total  annual  operating  energy  and  the  ethylene  glycol  velocity. 

The  total  annual  gross  recovered  increases  as  the  ethylene  glycol  ve- 
locity increases.  The  total  annual  gross  recovery  is  directly  proportional  to 
the  product  of  the  overall  heat  transfer  coefficient  and  the  total  recovery  unit 
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area.     The  total   recovery  unit  area  is  constant,  and  the  overall  heat  transfer 
coefficient  is  directly  proportional   to  a  fractional  power  of  the  ethylene 
glycol   velocity  (VEG      )  as  shown  in  (5),  hence  the  above  result  is  satisfied. 

The  total  annual   net  gross  recovered  increase  as  the  ethylene  glycol 
velocity  increases.     The  total   annual   net  gross  recovered  is  the  difference 
between  the  total  annual   gross  recovered  and  the  total  annual  operating 
energy. 

The  ethylene  glycol  velocity  equal  to  8  is  considered  as  the  maximum 
ethylene  velocity  for  the  system  to  run  under  low  pressure.  The  ethylene 
glycol  velocity  equal  to  8  is  taken  as  the  standard  value  for  the  system. 
Above  this  velocity,  the  system  is  considered  working  under  high  pressure 
which  requires  high  cost  fabrication  and  components. 

Figure  13  shows  that  the  total  annual   gross  recovered  increases  as 
the  number  of  rows  in  coil   depth  increases,  because  the  total   recovery 
unit  area  increases  the  number  of  rows  in  coil  depth  increases.     However, 
the  numerical  results  show  that  the  by-pass  was  open  for  preventation 
of  frost  formation.     This  limited  the  total  annual   gross  recovered,  as  noted 
in  the  figure,  by  reducing  the  ethylene  glycol  flow  through  the  make  up 
coil . 

The  total  annual  operating  energy  curve  shows  that  it  increases  as 
the  number  of  rows  in  coil   depth  increases,  because  the  total  annual  oper- 
ating energy  depends  on  the  air  pressure  drop  and  the  air  pressure  drop  depends 
on  the  number  of  rows  in  coil   depth. 

The  total  annual   net  gross  recovered  curve  shows  an  increase  until 
the  number  of  rows  in  coil   depth  equal  to  4,  then  decreases  as  the  number 
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of  rows  increases  beyond  4.  The  reason  is  that  the  total  annual  net  gross 
recovered  is  directly  proportional  to  the  total  annual  gross  recovered  and 
inversely  proportional  to  the  total  annual  operating  energy.  The  magnitude 
of  the  total  annual  operating  energy  gradually  increases  until  the  number 
of  rows  equal  to  4;  for  values  of  the  number  of  rows  greater  than  4,  the 
increase  is  much  steeper.  For  this  reason,  the  value  of  the  total  annual 
net  gross  recovered  tends  to  decrease  after  the  number  of  rows  equal  to  4. 

Figure  14  shows  the  effect  of  the  number  of  rows  in  coil  depth  vari- 
ation on  the  annual  net  energy  cost  recovered.  The  annual  net  energy  cost 
recovered  for  oil,  gas  and  electric  curves  increase  until  the  number  of  rows 
equal  to  4.  Further  increase  in  the  number  of  rows  results  in  fall  of  the 
above  values.  This  is  because  the  annual  net  energy  cost  recovered  is  directly 
proportional  to  the  total  energy  gross  recovered  cost  and  inversely  to  the 
total  operating  energy  cost,  as  already  explained  in  the  discussion  of 
figure  8.  The  significant  effect  of  the  total  operating  energy  cost  on  these 
curves  produces  this  result. 

Figure  15  shows  the  effect  of  the  number  of  rows  on  the  payback  period. 
Because  the  payback  is  inversely  proportional  to  the  capital  equipment  re- 
duction cost,  the  total  annual  net  energy  cost  recovered  and  directly 
proportional  to  the  total  initial  cost,  these  curves  increase  with  dif- 
ferent slopes. 

In  short,  figures  13  to  15  indicate  that  the  highest  annual  net  energy 
recovered  and  the  highest  annual  net  energy  cost  recovered  are  obtained 
for  the  number  of  rows  equal  to  4,  and  the  quickest  pay  back  period  results 
for  the  number  of  rows  equal  to  2.  Hence  the  value  of  the  number  of  the 
rows  equal  to  4  is  selected  as  standard  coil  dimension. 
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Figures  7  through  15  are  based  on  a  geographic  area  with  4000  to 
5000  degree  days. 

Figures  16  to  24  demonstrate  also  the  standard  values  of  the  fin 
spacing,  face  velocity,  ethylene  glycol  velocity  and  coil  depth  for  a 
geographic  area  with  5500  to  6500  degree  days.  This  set  of  conditions 
shows  the  same  results  regarding  the  four  parameters  as  in  the  first  set 
of  degree  days . 

The  only  difference  between  the  two  sets  is  that  the  annual  net  energy 
recovered  and  the  annual  net  energy  cost  recovered  have  greater  values  in 
the  second  set  than  in  the  first  set  and  the  pay  back  periods  have  lower 
values  in  the  second  set  than  in  the  first  set.  This  occurs  because  the 
total  annual  net  gross  recovered  is  directly  proportional  with  the  degree 
days.  The  same  influence  is  noted  on  the  annual  net  energy  cost  recovered 
for  oil,  gas  and  electric.  The  pay  back  periods  are  inversely  proportional 
to  the  annual  net  energy  cost  recovered  which  is  directly  proportional 
with  the  degree  days. 
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OUTDOOR    TEMPERATURE 


Figure  3.  Climatological  Data  for  an  Area  with  4000  to 

5000  Degree  Days,  The  outdoor  temperatures  are 
related  to  the  hours  of  occurrence  of  the  different 
months  of  the  heating  season.  After  U.S.  Dept.  of 
Commerce  Weather  Bureau,  Asheville,  N.C. 
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OUTDOOR  TEMP.- F 


Figure  4.  Climatological  Data  for  an  Area  with  4000  to  5000 

Degree  Days.  The  outdoor  temperatures  are  related  to 
the  hours  of  occurrence  of  the  different  months  of 
the  cooling  season.  After  U.S.  Dept.  of  Commerce 
Weather  Burea,  Asheville,  N.C. 
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OUTDOOR  TEMPERATURE 

Figure  5.  Climatological  Data  for  an  Area  with  5500  to  6500 
Degree  Days.  The  outdoor  temperatures  are  related 
to  the  hours  of  occurrence  for  the  different  months 
of  the  heating  season.  After  U.S.  Dept.  of  Commerce 
Weather  Bureau,  Asheville,  N.C. 
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OUTDOOR  TEMP.-F 

Figure  6.  Climatological  Data  for  an  Area  with  5500  to  6500 
Degree  Days.  The  outdoor  temperatures  are  related 
to  the  hours  of  occurrence  of  the  different  months 
of  the  cooling  season.  After  U.S.  DeDt.  of  Commerce 
Weather  Burea,  Asheville,  N.C. 


29 


OL  *  "id   'bS  "  V  'V3«V  UNO  AU3A033U  "1V101 
1  III''' 1 1 


r-~  lo 


r-         oi 


ro 


C\J 


•ha  /  9ol  *  ma  -  30vi  'akom3  OHiivaadO  ivioi 
i       i       iii'1' I 1 


•MA  /  „0L  x  fll8  -  M9NV1  P"?  M9V1  'A9M3N3  "IVIOi 


30 


"HA  /c 01  x  $  -  3UN3   '('0313)   AU3A0D3a  13 N 


I 


J 


LD  *3"  CO  CM  r—  O 

o        o       p        o        o        o 

CSJ  CM  C\J  C\J  CSJ  CM 


OS 

CD 

Is* 

10 

LO 

CI 

PI 

CTi 

cr> 

en 

•«A  /  ,01  X  $  -  a9N3  "(SV9)  AH3A033a  13N 
i         P        1         i  I         I         I I 1 

rocMr—  oo^cor^.  i£>m 


*HA  /  20L  x  $  -  0HN3   '(HO)   AcJ3A033a  13N 


1 

1 

1 

| 

| 

| 

SW3A  - 

|         1 

3HN3 

1 

'(•3313) 

L       ' 

XOVaAVd 

CO 

r-* 

lO 

s 

ro 

CNJ 

,— 

o 

CN 

* 

* 

ro 

ro 

ro 

ro 

en 

CO 

ro 

ro 

C\J 

cnj 

CNJ 

CM 

CXJ 

M 

CXJ 

CNJ 

CnJ 

CNJ 

CNJ 

CNJ 

CNJ 

1 

1 

| 

| 

| 

1 

,  i 

SHV3A 

1 

-  U9N3   * 
f          1 

(SV9) 

1 

xovaxvd 

CO 

r-. 

2 

ro 

CM 

o 

ON 

eg 

N. 

«3- 

sf 

** 

9 

«st 

•=r 

>» 

«sfr 

suv3A  -  ami  '(no)  »ovaAVd 


32 


m        oj 


-0L  *  •Id   -t)S  -  V  'V3bV  UNfl  AH3A003H  "1V101 


i 


av3A  /  goi  x  ma  -  sovi  'A9H3N3  9Niiva3dO  "ivioi 


1    1 

1 

1 

| 

! 

t 

1 

i 

..,!._ 

1 

| 

o 

C?i 

CO 

rs. 

UD 

LT) 

^J- 

CO 

CM 

o 

en 

CO 

r--. 

UD 

<jD 

lO 

kjD 

IT> 

U3 

to 

10 

lO 

UD 

" 

LT) 

HV3A  /  n0L  x  OLS  -  H9NVI  P"B  «9V1  'A9«3N3  "IVIOI 


53 


•MA  /  c0l  X  $  -  3UN3  'C3313)  AU3A003H  13 N 
■    i    '    I    '    '    I    I    I 


"MA  /  ,01  *  $  -   H9N3  *(SV9)  AH3AO03H  13N 


"UA  /  20L  x  S  -  0UN3 


'(110)  AH3A033M  13M 


34 


■  01  x  'Id  "OS  -  V  'V3W  UNO  AH3A003y  IVIOI 


■U  I    OL  *  flia  -  30V1  'A9H3N3  9NIlVM3dO  1V101 


•ma  /  ol  x  ma  -  hbnvi  Pue  asvi  * A9M3N3  hvioi 


35 


I L 


,co  r-»         <<o         ld  «3-  ro  cm         ■ —         o 

-01  x  '13   "OS  -  V  'V3HV  UNfl  AH3A003H  "1VJ.01 
_I 1 I ! I ! I I I 


OD 


■u  I  goi  x  ma  -  3ovi  'a9H3N3  9nuvh3<K)  ivioji 


1 

! 

1 

! 

I 

I 

| 

I 

1 

1 

CSJ 

o 

CO 

•x> 

"3- 

CNJ 

O 

CO 

ID 

<a- 

C\J 

o 

OS 

r-. 

r«s 

lO 

tc 

U3 

U3 

KO 

Lf) 

in 

irs 

» 

tO 

■a- 

•«A  /  „0L  x  nia  -  H9NV1  Pub  a9Vi  'A9U3N3  "IViOl 


56 


♦-.     «r-     Lfl 


r*. 

vo 

LT> 

^t-          ro          eg          r— 

o 

en  03  r^  ^o  tn 
o         o         o        o        o 

: 

■ 

**iA  /  r0L  x  $   - 
!            !            1 

3HN3 

'("3313)  AH3A033U  13N 
1 1 L__l 1 

"HH  I   Z0L  x  %  -   H9N3  '(SV9)  AH3A003H  13N 


'HA  /  20L  x  S  -  0HN3  '("110)  AH3A003*)  13N 


37 


SW3A  -  38AVd   '('3313)  >IOV9AVd 
J I I I I I I 


ooicor-.v£>un»fl-coe\Jr— 


SW3A  "  39AVd  '(•3313)  >13V8AVd 


oor^  id         lo         ^oocsjf—         ocrt 


sav3A  -  oaAVd  '(no)  >3vaAVd 


3S 


.01  *  "jj  'bs  -  v  'vaav  iiK'n  ambaoohm  "ivioi 


_L 


J 


CO       r— 


•ma  /  g0L  x  ma 


30V1  'A9H3N3  ONIlViBdO  1V101 
I I 1 I I 1 


•UA  /  „0l  *  ma  -  aSNVJ.  P"e  U0V1  'A9M3N3  1V101 


39 


....    ! 

1  1      1 

O 

t/i 

OJ 

n 

L 


> 


,_         _<  4->  O 


'HA  /  £0L  x  %  -  3UN3  'CD313)   AH3A033H  13N 


•»A  /  ,01  *  $  -H9N3   '(SV9)   AU3A(raH  I3N 
I        f         I         i         I         i         i         I         i 


I—  o 


•HA  /     01  ^  $   -  0HN3   '(110)   AU3A003H  13N 


40 


SHV3A  -  39AVd  '('33131  >IOV9AVd 
'         j         I         I I I 1 


1      1 

1 

SW3A 
I          1 

-  saxvd  ' 

SV9) 

>i3V3AVd 
I          i 

O             01 

cc 

IX)              LD 

*a-         ro 

F—                      O 

•X>            (X> 

sav3A  -  oaAVd  '(no)  >iovaAVd 


O  M  f— 


OL  *   'Id   ebS  -  V  6V3W  UNfl  A^3A003y   1V101 


CM  .— 


•u  I  g0L  ^  ma 


i 


30V1   'A9M3N3  9NIlVH3dO  "IVJ.01 
J 


J_ 


av3A  /  Ooi  x  nia  -  h9nvi  p"e  asvi  'A9H3N3  ivioi 


42 


—  or;  T. 


CO 

K 

<x> 

in 

«• 

ro 

C\J 

i- 

o 

CTi 

CD 

r*. 

r-*. 

p-. 

r-^ 

r»* 

r^ 

r»» 

(•*» 

r^ 

r^ 

VD 

10 

u? 

1 

! 

| 

aA  / 

1 

e0L 

1 

X 

s  - 

1 

aaN3 

'C 

! 

3333) 
1 

Ab3A003a 

1 

JJN 

tN, 

to 

LO 

^ 

ro 

CM 

o 

o 

o 

o 

o 

•MA  /  70l  x  $  -  b9N3  '(SV9)  AH3A003a  I3N 

I       3       I       I       I I I I— J 

<a-nc\j         r—         o        cti         co         r-.        lo 


•aA  /  20L  x  S  -  0aN3  '(310)   Aa3A003a  13N 


43 


CM      r- 


*HA  /  „0L  x  0L9  -  H9NV1  Pue  H9V1  'A9U3N3  IVlOi 


44 


OJ  r— 


„0L  x  -1J  "OS  -  V  'V3bV  UNfl  AH3A0D3b  TVlOi 
I         i         i         i         I         I         I         I I 


t-t       *»-       in 


'HA  /  g0l  x  ma  -  30V1  'A9U3N3  9IJIlVa3dO  "1V101 
l         I         I         I 1 1 I I !___! 


•«A  /  n0L  x  Ilia  -  H9NV1  PUE  «9Vi  'A9H3N3  "IV101 


45 


00 

1 

| 

CO 

1 

| 

r —        r —         r — 

■UK  /  r0L  x  $ 
1      bl         ! 

r-*            VD             UD            (X)            ^O            U3 

-  3bM3  '(3313)   Ab3M)03a  13N 
i         <         1         1         1         1 

tn 

m 

ir> 

IX)            ix)          IX) 

IX)             LT)            LO             LO            tO             1X1 

'HA  /  z0l   x  $  -  H9N3  '(SV9)  Ab3A003b  i3N 
i    I    I    I    i    I    I    I 1 


'HA  /  ,01  *  $  "  0bN3  '(110)  Ab3A003ci  13N 


45 


1        1 

CO 

1 

1 

1 

| 

CO            CN) 

SW3A  - 

1      1 

CM 

38AVd 

I 

cm         cm         cm         r- 

T3313)  >l3Va  AVd 
1         1 J 1 

■a-        cm 

© 

CO 

Ifi 

<tf 

CM          O 

00 

t£)             «*            C\)            O 

SUV3A  -  98AVd  '(SV9)  'MDV9  AVd 
I    I    I    'I    I I 


■—  r—  O 


ScW3A  -  OSAVd   '(110)   >UV9  AVd 


CHAPTER  5 

"HEAT  PIPE  RECOVERY  UNIT" 
General   Description 

figure  25     shows  the  typical   ventilation  system  using  the  heat  pipe 
recovery  unit. 

The  efficiency  of  a  "neat  pip?"  thermal   recovery  ur.it  is  The  measure 

of  its  ability  to  recover  thermal   er-ergy  from  one  air  stream  and  to  transfer 


Exhaust- 


Supply 


..  1 


Exhaust  Fan 


•  2 

^ 

•  3 

& 

Heat  Recovery 
Unit 


Supply  Fan 


Figure  25.  Arrangement  of  Heat  Pipe  Recovery  Unit. 

it  into  another  air  stream  flowing  in  the  opposite  direction.  The 

recovery  factor  (R)  is  the  percentage  of  availanie  sensible  thermal  energy 
that  is  recovered.  The  recovery  factors  are  based  on  sensible  neat  transfer 
only,  even  though  latent  heat  is  transferred  when  condensation  occurs  on  the 
unit. 
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The  required  information  for  analysis  is  supplied  by  the  manufacturer 
of  the  unit.     The  data  used  here  is  from  Q-DOT,  but  it  is  typical  of  other 
manufacturers  of  heat  pipe  units  and  even  of  data  for  other  types  of  exchangers 
(heat  wheel,  recuperator,  etc.) 

The  recovery  factor  of  a  unit  varies  with  the  geometry  of  the  unit. 
For  the  heat  pipe  type  the  variables  are  the  number  of  fins  per  inch,  the 
number  of  rows  of  heat  pipes,  the  mass  flow  ratio  of  the  two  streams,  and  the 
air  density. 

The  pressure  drop  across  each  side  of  a  thermal  recovery  unit  is  de- 
pendent upon  the  number  of  fins  per  inch,  the  number  of  rows  of  heat  pipe, 
the  air  velocity,  and  the  air  density. 

Analysis  can  be  made  by  using  characteristic  curves.     For  this  analysis 
characteristic  curves  for  Q-PIPE  thermal   recovery  units  produced  by  Q-DOT 
CORPORATION  [13]  were  chosen. 

Air  density  corrections  due  to  temperature  and  moisture  variations  are 
negligible.     Altitude  density  corrections  are  significant  and  are  included. 
Temperature  and  moisture  density  corrections  must  be  considered  for  appli- 
cations other  than  those  within  the  comfort  range. 

No  control   is  required  for  summer-winter  switch  over.     Because  the 
heat  pipe  units  are  bidirectional,  transferring  heat  from  the  high  temper- 
ature side  to  the  other  side,  they  automatically  switch  to  the  precooling 
mode  whenever  the  exhaust  becomes  cooler  than  the  air  in  the  supply  duct. 

Theory  and  Modeling  Equations 
Heat  Balance  and  Transfer  Equations 

The  heat  transfer  and  balance  equations  for  this  system  were  written 
as    (see  figure  25): 


VT1  +  [Ra(YV] 

VV[VYT1)] 
T3«T1+[Rb(T2-T1)] 


T4*Tr^VTi» 
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(27) 

(28) 
(29) 
(28a) 
(29a) 


Equations  28  and  29  are  used  when  supply  air  flow  is  smaller  than 

exhaust  flow.     Equations   28a  and  29a  are  used  when  supply  air  flow  is 
larger  than  exhaust  flow.      If  air  flows  are  equal    {P.&  '  Rb),  any  of    these 
equations  may  be  used. 

The  theory  cresen+ed  "'s  that  developed  by  Isotherrnlcs,   Inc.   [14]. 
Details  of  determination  of  R    and  R.    are  given  in  Appendix  C. 

Frost  Control  Determination 

The  neat  pine  units  are  easily  adapted  to  face-and-by  pass  control, 
for  frost  prevention  and  temperature  control,  as  shown  in  figure  26. 


Exhaust 


Supply 


Heat  Recovery  Unit 

Exhaust  Fan 


Supply  Fan 


Figure  25.     Schematic  Diagran  for  Frost  Control, 
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The  preceding  schematic  diagram  is  followed  to  prevent  frost  in  the 
exhaust  air  side.     If  the  exhaust  air  temperature  T4  <  32     F,  part  of  air 
supply  is  bypassed;  this  option  is  built  in  the  program. 

Another  possibility  to  prevent  frost  is  to  pre-heat  the  supply 
air  before  it  enters  the  recovery  unit.     This  was  not  considered  here 
because  of  its  high  initial  and  operation  cost. 

Characteristics  Curves 

A  6  row  unit,  14  fins  per  inch,  is  chosen  for  this  analysis.     The  per- 
formance curves  are  shown  in  figures  27,  28,and  29,  and  table  1.     The  general 
form  of  these  curves    was    found  as: 

VS  =  [-7.3041216  +  1.027933  AV  -  0.253238  x  10"4  AVZ] 

+  [0.5277119  x  10~2  -  0.6327463  x  10~4  AV  +  0.280454  x  10"7  AV2] 

H  +  [-0.6582271  x  10"6  +  0.3537444  x  10"8  AV  -  0.3220021 

x  10"11   AV2]  H2  (30) 


R  ■  [0.4819098  +  0.3368778  G  -  0.6189942  x  10"1  G2] 

+  [-0.366823  x  10"3  +  0.1159242  x  10"3  G  -  0.87691  x  10"4  G2] 

VS  +  [0.2868927  x  10"6  -  0.2522431  x  10~6  G  +  0.1094314  x  10"6 

G2]  VS2  (3D 
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APD  =  [-0.7787704  x  10"1   -  0.1196543  x  10"4  H  +  0.1090248 
x  10"8  H2]  +  [0.8309141  x  10"3  +  0.9480123  x  10"7  H  - 
0.6518396  x  10"11   H2]  VS  +  [0.1757011  x  10"5  - 
0.2584665  x  10"10  H  +  0.1227925  x  10-13  H2]  VS2  (32) 


F     ■  1.08  -  3.825  x  10"5  H  +  3.75  x  10"10  H2  (33) 

Details  of  equations  30,  31,  32, and  33  are  given  in  Appendix  C. 


Energy  Equations 

Energy  savings  at  design  temperature  may  be  calculated  as  follows, 
assuming  100  percent  energy  utilization  efficiency: 

q  =  1.08  cfms  x  F,.  x  (Tg-T.,)  (34) 

Regarding  the  energy  recovered  per  season,  data  from  figures  3,  4,  5, 
and  6  have  been  used  in  the  calculational  program. 

The  fan  operating  energy  is  as  in  equation  11,  and  the  net  energy 
recovered  is 

q  =  q  (energy  recovered)  -  q  (operating  energy  for  fans)      (35) 

The  Cost  Equations  -  Heat  Pipe  Units 

To  determine: 

The  Total  Initial  Cost: 

cost  ($)  =  I   cost  (Heat  pipe  unit  +  Labor  +  Fans)  (36) 


s: 


The  Fan  Operating  Energy  Cost: 

cost  ($)  =  (Fan  BHP)  x  0.746  x  ^^  (Elec.  cost)  x  y—  (16) 

The  Total   Capital   Equipment  Reduction 

Eguations  19,  20  and  21  are  also  used. 
The  Total   Energy  Gross  Recovery: 

Equations  22,  23  and  24  are  also  used. 
The  Net  Energy  Recovery  Cost: 

Equation  25  is  used. 

The  Years  of  Pay  Back: 

Equation  26  is  used. 

The  initial   cost  was  selected  on  the  basis  of  existing  data,  and  may 
be  changed  for  more  recent  information. 

Calculational  Scheme 

A  computer  program  (shown  in  Appendix  D)  was  written  for  finding  the 
equations  of  each  curve  in  figures   (27),    (28),  and  (29).     By  selecting  any 
three  arbitrary  points  on  these  curves,  R  «  f(VS),  VS  =  f(H),  and  APD  =  f(VS) 
are  obtained.     The  same  program  was  then  used  to  find  the  general   equation  of 
these  curves  as  function  of  both  variables  for  each  figure   (i  .e.,  equations 
30,  31,  and  32). 

A  computer  program  (shown  in  Appendix  E)  was  written  which  by  inputting 
equations  30,  31,  and  32,  air  flow  rates,  air  temperatures,  fan  efficiencies, 


53 


system's  component  cost.and  unit  heating  and  cooling  costs,  allowed  cal- 
culation of  R    and  R.  .     With  R     and  R.    known,  the  program  then  determined 

8  *D  3D 

T3  and  T4  as  in  equations  (28)  and  (29)  or  (28-a)  and  (29-a).  The  program, 
using  the  selected  input  values,  determines  if  the  resulting  T"4  is  below 

32  °F. 

If  T.  is  below  32  °F,  the  program  option  is  selected  which  allows 
part  of  the  supply  air  (as  shown  in  figure  26)  to  be  by-passed;  then  by 
fixing  T.  to  be  equal   to  32  °F,  T3  and  then  T"5  were  found. 

The  program  then  determines  the  energy  recovered  per  season,  the 
operating  energy,  the  cost  results, and  the  pay  back  period. 
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STANDARD  AIR  VELOCITY,  VS  -  FPM 
Figure  27.  Recovery  Performance  for  6  Row  Heat  Pipe  Recovery  Unit, 
14  Fins  Per  Inch. 


w*m. 


,m^r^mm& 
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u         STANDARD  AIR  VELOCITY,  VS  -  FPH 
Figure  28.  Elevation  Correction  for  6  Row  Heat  Pipe  Recovery  Unit, 
14  Fins  Per  Inch. 
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STANDARD  AIR  VELOCITY,  VS  -  FPM 
Figure  29.  Air  Friction  for  6  Row  Heat  Pipe  Recovery  Unit, 
14  Fins  Per  Inc. 


ELEVATION 

F« 

ELEVATION 

F« 

Sea  Level 
2000  It. 
4000  It 

1.080 

1.005 
0.934 

5000  ft. 
6000  ft 
8000  ft. 

0899 
0.863 
0.798 

Table  1.  Air  Flow  Correction  Factor. 


CHAPTER  6 

SPECIFIC  SYSTEM  ANALYSIS 
Results  and  Discussion 

For  the  heat  pipe  system,  the  effect  of  air  flow  velocities  on  the 
system  performance  has  been  studied. 

Figure  30,  which  based  on  an  geographic  location  with  4000  to  5000 
degree  days,  shows  the  effect  of  variation  on  the  air  face  velocity  (V) 
on  the  total  annual  gross  recovered  (TAGR),  the  total  annual  net  gross  re- 
covered (TANGR),  the  total  annual  operating  energy  consummed  for  the  fans 
(TA0E)  and  the  face  area  of  the  recovery  unit. 

The  total  annual  gross  recovered  decreases  as  the  air  face  velocity 
increases,  because  the  total  annual  gross  recovered  depends  on  the  inlet  and 
outlet  air  fresh  temperature  T1  and  T3  correspondingly,  provided  T1  is  given 
(constant).  T.  is  variable  and  depends  upon  recovery  factor  (R).  The 
recovery  factor  decreases  as  the  air  face  velocity  increases,  because  of 
smaller  area  (see  figure  27;  recovery  performance)  for  constant  air  flow 
capacity. 

The  total  annual  operating  energy  curve  shows  an  increase  as  the  air 
face  velocity  increases,  due  to  air  pressure  drop  (APD). 

The  total  annual  net  energy  recovered  curve  decreases  as  the  air  face 
velocity  increases,  because  it  depends  upon  the  total  annual  energy  recovered 


sc 
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and  inversely  dependent  on  the  total  annual  operating  energy. 

Recovery  unit  face  area  curve  shows  a  decrease  as  the  air  face  velocity 
increases,  if  the  air  flow  capacity  is  constant. 

The  economical  range  of  velocities  is  in  between  400  to  600  fpm. 

The  total  annual  gross  recovered,  the  total  annual  net  gross  recovered, 
and  the  recovery  unit  face  area  are  higher  at  low  air  face  velocity  which 
means  more  expensive  recovery  unit  but  the  total  annual  operating  energy 
is  low;  however  the  reverse  is  true  at  higher  velocity. 

Figure  31,  which  based  on  an  geographic  location  with  5500  to  6500 
degree  days,  give  the  same  result  and  conclusion. 

The  main  difference  between  set  1  and  set  2,  degree  days  is  that 
the  total  annual  gross  recovered  and  the  total  annual  net  gross  recovered 
are  directly  proportional  to  the  degree  days. 
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CHAPTER  7 
CONCLUSIONS  AND  RECOMMENDATIONS 

The  purpose  of  this  study  was  to  analyze  two  heat  recovery  units,  in 
order  to  determine  the  influence  of  specific  parameters  on  the  energy  re- 
covered.    Thereafter,  it  was  to  decide  the  best  system  suited  for  a  specific 
purpose  in  order  to  conserve  some  energy  being  wasted. 

The  basic  equations  describing  the  two  systems,  run-around  coil  system 
and  heat-pipe  recovery  unit,  were  developed  and  the  economic  considerations 
were  investigated. 

The  reason  for  selecting  these  two  types  was  to  provide  analysis  for 
two  basically  different  systems  as  shown  below, 

A.  Run-Around  Coil  System: 

1.  It  is  indirect  heat  transfer  system. 

2.  Several  coil  parameters  were  varied  and  subjected  to  the 
study.     Those  coil  parameters  which  cause  the  variation 
of  heat  transfer  area. 

B.  Heat  Pipe  Recovery  Unit: 

1.  It  is  direct  heat  transfer  system. 

2.  The  characteristic  curves  of  this  type,  which  were  submitted 
by  the  manufacturer,  was  subject  for  the  study  (no  parameters 
were  selected,  i.e.  fixed  coil  dimensions). 

3.  According  to  the  study  of  this  type  same  can  be  applied  for 
Rotary  and  Plate  recovery  types  using  characteristic  curves 
from  the  manufacturer. 
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The  study  for  run-around  coil  system  was  done  by  "built  up"  system  of 
standard  components  available  from  a  variety  of  manufacturers  and  notice 
the  effect  of  variation  of  each  parameter,  keeping  the  others  constant, 
on  the  net  energy  recovered  (savings)  and  the  coil  size  (initial  cost)  and 
repeat  it  for  all  the  parameters. 

The  analysis  for  heat-pipe  recovery  unit  depends  on  performance  data 
reported  by  manufacturer.     The  effect  of  air  flow  velocity  variation  on  the 
net  energy  recovered  (saving)  and  the  coil  size  (initial   cost)  was  studied. 

The  study  for  the  run-around  coil  system  shows  that: 

1.  Increasing  the  number  of  coil   tubes  or  decrease  the  fin  spacing 
leads  to  increase  in  coil   cost  and  the  pressure  drop.     It  increases  the 
heat  transfer  area  (increase  in  initial   coil   cost).     For  this  reason  energy 
recovered  increases  but  the  net  energy  recovered  (saving)  shows  maximum 
value  for  a  specified  values  of  the  parameter. 

2.  Increasing  the  air  flow  velocity  decreases  the  energy  recovered 
and  increases  the  pressure  drop.     The  reduced  coil  size  needed  for  the 
specified  air  flow  quantity  reduces  the  net  energy  recovered. 

3.  Increase  in  the  ethylene  glycol  velocity  increases  the  heat  transfer 
coefficient.     This  increases  the  energy  recovered,  and  increase  glycol 
pumping  cost  (operating  energy  cost)  which  shows  that  effect  on  the  net 
energy  recovered  (saving).     A  noticeable  increase  in  net  energy  recovered 
(saving)  was  shown  but  glycol   velocity  is  limited  to  8  fps  for  the  system 

to  be  operative  under  low  pressure. 
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The  study  for  heat-pipe  recovery  shows  that  the  energy  recovered  de- 
creases as  the  air  flow  velocity  increases.     A  decrease  in  face  coil  area 
causes  an  increase  in  the  operating  energy.    Thus  a  decrease  in  net  energy 
recovered  was  noticed. 

The  problem  can  be  extended  in  different  directions.     The  effect  of 
latent  heat  was  not  considered  because  the  sensible  heat-gain  is  much 
greater  than  the  latent  heat  gain  for  most  applications.     The  initial  cost 
of  the  systems  and  their  accessories,  the  heating  and  cooling  cost  were 
selected  to  be  reasonable  and  maybe  changed  to  reflect  further  increase. 
The  additional  cost  necessary  for  pump  and  fan  size  increases  was  not  con- 
sidered here.     This  should  be  added  in  a  more  detailed  analysis.     No 
maintenance  cost  is  included.     It  varies  considerably  with  the  unit  size, 
proficiency  of  the  servicing  organization  and  the  operating  hours.     Mainten- 
ance cost  should  be  included  if  it  is  available. 

The  economic  feasibility  of  inserting  a  unit  in  an  existing  system 
depends  primarily  on  the  savings  in  operating  costs.     These  savings  are 
dependent  on  the  required  volume  and  temperature  delivered  to  the  system, 
the  volume  and  temperature  of  the  exhaust,  hours  of  operation,  cost  of 
energy,  and  the  local  design  temperature. 

From  these  calculations,  it  can  be  seen  that  either  of  the  two  heat 
recovery  systems  can  be  installed  with  a  little  or  no  extra  cost.     If 
savings  in  heating  and  cooling  plants  costs  are  taken  into  account,  con- 
siderable savings  could  be  availed  in  operating  costs. 
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Appendix  A 
Derivation  of  Heat  Balance  and  Transfer  Equations 
The  terms  which  are  used  in  equation  (2)  are 
Q,  =4.5  cfms 

Q~  =  4.5  cfme 

G,  -  C   ■  Qi  =  0.24  x  4.5  cfms  =  1.08  cfms 

1  pa   1 

G„  =  C   •  Qo  *  0.24  x  4.5  cfme  =  1  .08  cfme 

2  pa   2 

From  equation  2 

S  <VV  ■  Gi  <W  (A_1) 

(VT,)  -  (T6-T,) 

=  U1A1    ,n[(T5-T3)/(VTl^         '"  C0Unt6r  fl°W 

(A-2) 

Solving  for  T,  in  equation  (A-l),  and  substituting  into  equation 
(A-2).  The  corresponding  value  of  Tg  is 


Rt^-e 
where 

D,   =  U1A1    [1/G3  -  1/S,] 

and 


(A-3) 


RM1   -  S,/G1 


Si 
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If  Tc  is  an  estimated  value  to  start  the  iteration,   its  new  value 

0 


after  finding  Tg  in  equation  (A-3)  is, 


RM2-e 
where 

D2  '  U2A2  [1/G4  -  1/G2] 

and 


RM2  =  G4/G2 

Equation   (A-4)   is  used  to  iterate  on  TV  if  the  by-pass  value  is  not 
open.     If  not, 

D2 

RM2-e  c 

Iteration  is  continued  until  a  temperature  balance  is  reached,  so  that 
from  equation  (2), 

VV  <VV  mv  (A"6) 

T4  =  T2  -  (T5-T6)  RM2   if  by  pass  value  is  closed,      (A-7) 


and 


T4  =  T2  "  'T5"T7^  RM2   if  by   pass  value  is  0n'         ^A"8^ 
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Frost  Control   Determination 

The  program  runs  the  selection  and  checks   to  see  if  the  resulting 
Ty  is  below  32  °F  as  follows: 

Since  the  study   is  to  find  a  maximum  heat  recovery  as  in 
equation   (A-i),  by  analyzing  the  left  term  it  is  found  that  T5  is  constant. 
For  a  maximum  heat  recovery  Tg  is  selected  to  be  equal   tc  T-j ,  as  initial 
value  for  the  iteration,   in  order  to  keep  a  maximum  value  for  trie  left 
term  of  equation   (A-l),  and  start  with  T_  »  32     F. 

Development  of  Hake  up  and  Exhaust  Air  Coils   Equation 
The  circuiting  coil   is  chosen  of  vertical  headers  feeding  horizontal 
tube  circuits   in  parallel,   as  shown  in  figure  A-l.     The  flew  of  ethylene 
glycol    through  the  U  bends    is  counter  to  tne  flew  of  air.     Pure 
counter  flow  was  assumed  in  the  coils  DI 

OD 


Section  A 
Figure  A-l.     Schematic  Diagram  for  Run-Around  Coil  toncem 
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Further  assumptions  are 


OD  -  I  "  «  0.052  ft. 


TK  =  0.017"  =  0.0014  ft. 


Sh  =  1  |  "  =  0.135  ft. 
S.  =  1  f  "  =  0.135  ft. 


OFD  =  2.44  x  OD  =  0.127  ft. 
FTK  =  0.014  "  =  0.0011   ft. 

The  values  of  OD,  TK,  Su  and  S£  were  selected  from  [5]  consistent 
with  equations  5  and  7. 

The  values  of  OFD  and  FTK  were  selected  to  maintain  clearance  of  0.1 
inch  between  every  two  adjacent  fins  [8]. 

One  of  the  most  important  steps  in  the  study  of  heat  transfer  is  to 
estimate  the  pressure  drop  and  the  rate  of  heat  transfer  for  turbulent 
flow  across  tube  bundles.     But  before  such  an  estimation  can  be  made,  the 
Reynolds  number  must  be  determined.     There  are  cases  where  laminar  or 
transition  flow  takes  place  in  the  passage  between  tubes.     In  most  actual 
applications,  however,  the  flow  is  turbulent. 

The  tubes  may  be  arranged  either  in  staggered  form  or  in  in-line 
form;  only  the  latter  one  is  considered  here.     To  achieve  more  generality, 
it  is  considered  that  there  is  no  change  in  the  tube  diameter  or  the 
spacing  from  one  bundle  to  another. 

From  equation  (6),   the  pressure  drop  in  cross  flow  was 


70 


f  G2  L  D     °'4     S     °-6 

APD  =  — K£ (**■)         (f1)         lb/in2  (A-9) 

(5.22xlOIU)   DyS  *s       sh  sh 

developed  by  Gunter  and  Shaw  [9]  and  presented  by  Kern  and  Kraus  [10].     The 
numerical   factor  (27.69)  in  equation  (6)  is  to  convert  the  dimensions  to 
inches  of  water.     The  main  terms  in  equation  (A-9)  can  be  found  by 

n    .  4  FNV 
uv  "      FRS 


where 


FNV  =  [S,  x  HT  x  Z  -  |  0D2  x  Z  x  RN  x  W  -  |  (OFD2  -  0D2)  x  FTK  x  RN  x  W] 


FRS  =   (Sf+S  )   RN  x  Z 
Sf  ■  6n(0FD2-0D2)  FN 


SQ  =  12n  (1-  FN  x  FTK)  OD 


Z  =  CFM  /  (HT  x  VA) 


G     =£ 
s       as 


where 


as  =  [HT  -  {OD  -   (OFD-OD)  FTK  x  FN}   RN]  Z 


.  p     _  0.635978195 

"  62-4  Ta 


where 

.  14.7  x  144 

P  K-5     ^    v    T 


1      =  Sr   x  W 
P         £ 


(W    /     MS) 


0." 


*     can  be  considered  to  oe  1 ,   because  the  viscosity  variation 

between  the  caloric   temperature  and  the  tube-wall    temperature  is  small 

enough  in  the  application  considered  here  tnat  if  can  be  neglected. 

To  find  the  friction  factor  for  air  side   (f ) ,    it  is  necessary  to 

evaluate 

pV     D         3D 
^     .      ra    v  .     s     v 


Figure  A-2,   shows  correlations  presented  by  Kern  and  Kraus   LlOj. 
The  equation  of  the  line   is 


f  =  0.0025288  -  4.115  x 


in-9 


5     S  !C!  2       i    «       6 

Figure  A-2.     Transverse-Fin  Heat  Transfer  and  Pressure 
Drop,    (a)  Jameson   [Trans.   ASME,  67:623- 
642   (1945)];    (b)   Gunter  and  Shaw  [Trans.  ASHE, 
67:    c43   ;i9i5}]. 
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Equation  (7),  the  ethylene  glycol  pressure  drop  was 

1    75 
PEG  ■  1.495  (^p)    "      [0.0039  WZ  +  0.0875  (W-l)  +  0.3]  (A-10) 

as  given  by  Stoecker  [5].     Stoecker  [5]   recommends  addition  of  5  ft.  of  water 
head  for  interconnecting  piping,  when  a  glycol  velocity  varies  from  1  to  6  fps 
Since  the  head  loss  for  steady,   incompressible,  constant-area  flow  is 

h     ,  ,  ^  VEGi  =  K  VEGi  (A_ll) 

\       f  DI  2g  K  2g 

by  considering  3  fps  glycol  velocity  in  equation  (A-ll),   K  becomes 
equal   to  35.78  and  the  head  loss  in  equation   (7)   is  equivalent  to  5  ft 
water. 

Derivation  of  Ethylene  Glycol  Equations 
If  the  ethylene  glycol  with  following  physical   properties  (specific 
heat  and  specific  gravity) 
C     •  0.602 

S  *  1.1 

is  mixed  with  water  by  X  percent  by  weight,  so  the  mixture  physical 

properties  is 

C    «  (1  -  0.398  X) 

S  =   (1  +  0.1   X) 
ps  «  62.4  (1   +  0.1   X) 
then  these  relationships  lead  directly  to  equations   (8)  and  (9). 
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Derivations  of  Operating  Energy 


m   62.3  x  CFM  x  APD 
Fan  BHF  -  12  x  33000  -  £ff# 


0.000157  x  CFM  x  APD  (A-12) 

EfT. 


Fan  operating  energy 

=  (fan  BHP)  x  2545 


D    dud   GPM  x  S  x  8.333  x  PEG 
PumP  BHP  =  33000  x  Eff . 


GPH  x  S  x  PEG  (A-13) 

3960  x  Eff. 


Pump  operating  energy 

■  (Pump  BHP)  x  2545 
Equation  (A-12)  is  developed  by  [11]  and  equation  (A-13)  is  developed 
by  [12]. 

Derivations  of  the  Cost  Equations 
Equation  (14)  is  based  on  weight  of  coil  materials  as  presented  by 
Edward  [8].  A  multiple  of  this  cost  is  necessary  to  accomodate  material 
cost  increases.  A  factor  of  8  greater  than  Edwards  value  was  used  for  the 
results  presented  here.  However,  the  multiple  may  be  changed  to  reflect 
further  increases  or  more  accurate  cost  information. 
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Initial  costs  for  pump,  labor,  ethylene,  glycol  and  fans  in  equation 
(15)  were  selected  to  be  reasonable  and  may  be  changed  to  reflect  further 
increase. 

The  cost  recovery  and  operation  cost  of  the  system,  as  in  equations 
from  16  to  26,  were  based  on  the  following: 

a)  Heating  costs 

It  is  assumed  that  as  in  [7]  that 

oil   cost  ■  0.3      S/gallon  and  its  heating  eff.  ■  0.82 

gas  cost  =  1.05     $/1000  ft3  and  its  heating  eff.  ■  0.8 

electric  cost  =  0.025  $/Kw-hr  and  its  heating  eff.  =  1.0 
Then 

oil  heating  cost  =  0.26  $/Therm 

gas  heating  cost  =0.13  $/Therm 

electric  heating  cost  =  0.73  $/Therm 

b)  Cooling  Costs 

Refrigeration  system,   (i.e.  chiller)  so 

Electric  cooling  costs  =  0.025  Kw_hr  x  1     jQn 

"  °-025  ToT 

c)  Working  hours 

Working  hours  values  were  taken  from  figures  (3)  and  (4)  for  4000  to 
5000  degree  days  and  from  figures  (5)  and  (6)  for  5500  to  6500  degree  days. 

d)  Average  equipment  costs 

At  design  outside  temperature  (boiler  for  winter,  chiller  for  summer) 

the  average  equipment  cost  is  considered  125  =rr —  for  boiler,  and 

Thern 

90  =4-  for  chiller. 
Ton 
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The  values  in  sections  a,  b  and  d,  above  are  as  given  by  Trane 
Bulletin  [7],  but  may  be  changed  for  nore  accurate  cost  information. 
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fFMS=    SUPPLY    AIR    FLnjj     IN    CFJJ 
J'luf.     CYHAOST     MR     flHU        IN     CFM  _     _.., 

C    UT-    WINTFR    TF"pfraTIIPF     IN    nFRRFF    FAH. 
KtI    SIIMMFR    TEMPERATURE     IN    nFr,REE    FAH 


£  *      „«:   Si^PAtTNfi    WF'Sp.v'r.nsT 'for    the   FAN   IN    m 

,?  J      nip:    no  FRIT  I  MR    ENERGY    COST    FOR    THE    POMP     IN    ($) 

•  "l^l  6NERRY  SRflSS  p^nv"?T'N!rR,T^n  ,N  |HRS  /YR  ) 
JFFI  =  WINTFR  EQUIVALENT  HEATING  I  DAD  IN  '  "R|»'.L5» '. 
ilFl"    SIMMFR    EQUIVALENT    CnOLIfB    LOAD     IN     IHRS./YR.I 

nnf*  nn.  heating,  cost  in  is/tmwhi 


c  *********       _(  (   RIIN_tRn,,Nn  r.nn.  system     ))- 

CC  SC 

cc     :   ST 

£  *    vf&  ethylene  RiVdni,  vfj  nriTY  in  FPS. 

9-  *         VA«AIR    SinE    FACE    VEI.nr.ITY     IN    FPH. 

>  i   nn=  ouT-sinE  cnn.  I»«^  diamstm  IN  feet. 

r  -     FN-  mil   FIN  SPACTNr.  IN   FINS/INCH  . 

rc         *      kill  nit .  n<*  rows  in  hfight 

rC  :      PFrW-=   ?TnHY^NFRnR.SYr!S..npFlFSS..RF   HROP    IN    .FT.   nF   H20    ) 

>  *  FFi     FAN     FFFICIFNCY     IN      IS) 
£                  *  ■         PF=     PUMP     EFFir.IEwr.Y_INo(.I) 

c 

c 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

|  :  FT?rR,c:  ffitt'NtfHrfensnVnjt 

C  ;  |KF.  Sans   ^T^cnsy.}^ '  tion   „  ,„ 

\  I    Tfn£:    M    AN^^'-n^^T^.'^ER.Y.IN^RTU/YR.I, 

c 
c 
c 
c 
c 
c 

C  IMPLlcTT"  REAt.*R(  A-H.O-Z) 

SlM'FNSI,nNlnrnT(50l,ni.nHl50I.RNEWT(50).RNEHHI50) 
REAL    WTll/n.O/.STll/<)5.0/ 

REAn(5,4M    KK 
REARI5,*) 


liCP^ICPlUlCL.EICr,.EICF.DOT.Dr,T.DET.DPKW. 
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•  DPTR  ,nCHT.riCCT,V!Sl,VIS2 

5  RFA0(5.*)     VEG.VA1 ,VA2.FN1,FN2,W1.W2 
WT1-WT11 

IFIMM.FO.l  .AN0.MM.I.T.2)  r,n  TO  R 

IF(MM.c o.?.  AN0.MM.LT.3)  no  TO  9 

IF  (MX.F0.3.AN0.MM.I  T.4  I  r,0  TO  10 

IF(MM.F0.4.ANn. MM.LT.5)  GO  TO  11 

IFIMM.FO.5.AN0.WM.I.T.A)  GO  TO  13 

IF(MM.F0.6.  AND. MX. LT. 7)  GO  TO  46 

IF  (MX. to. 7  .ANn.MM  .1  T.fl  )  GO  TO  48 

IFfMM.FO.R. AND.MM.IT.9)  00  TO  53 

..  F  1  I  T  ^  (  L.        ~>  Q    \ 


R  WR1TFIA.2R 

GO  TO  3  5 
9  WRITF16.29) 

GO  TO  3  5 

10  WRITFI6.34) 
GO  TO  35 

11  WRITF(6,501 
GO  TO  35 

13    WRITF(6,5l) 

GO    TO    35 
46    WRITF|h,49) 

GO    TO    35 
4B    WRITF[fc,52) 

GO    TO    35 
53    WRITF(6,54) 
35    PRINT    19 

PRINT    22 

PRINT    19 

WR1TEI6.299I 
C 
0 

c 

011=001-2. *TKl 

012=002-2. *TK2 

0F01=2 .44*001 

0F02=2. 44*002 

HT1 =SH1*RN1 

HT2=SH?*RN2 

71=CFMS/(HT1*VA11 

72=CFMF/(HT2*VA2) 

W01=SI.1*WI 

wn2=SL2*W2 

AO  1  =  11  1  .ton  1*7  1  *RM*W1  1/(42.0) 


CFMS.r.  =  MF.WTl.WT2.STl.ST2.VFG.X.nni.002.VAl,VA2. 
FN1  .FN2.SH1  .SH? .SI  I .SI  2.RNI .RN2.W1 .W2.FE1.FE2.PE 
WF  =  l   ,SFF|   .MCR-^ICH.flCL.FICr-.FICF.TKl.T^.OO.. 
nr,T;oeT,OPKW,nRTR,OCHT,OCCT.VISl,VIS2,FTKl.FTK2 
l 


Alll  =  lll.*-'JilI~'_l-^','l""l''>—  c.wi 

A02=( 11 .*On2*7?*RN2*W2> / (42 .0 ) 

=  RS1  =  I  I '  132.0*7 ftRNl  1/(7.0)  )*IFN1*MOF01**2)-(OD1**2)> 

1FRsf'ul3?lot7?tRN?  f/17.01  )*IFN2*(  (  OF02**2  1-1  002**2  )  ) 
1    +2.*002*"(  l.-FN2*FTK2)  ) 


AI=FRS1*W1 

A2=FRS2*W2 

01=4.5*CFXS 

0?=4.5*CFMF 

G1=1.0«*C*MS 

G2=1.0R*CFME 

TAHGR=0.0 

TACGR=0.0 

NSET=l 

U.-1 


79 


q97    RF4D(5.1  )     INDEX 
REiO     (5,11    N 

12  ??aii2(  l"A  )' Vot  ( i  i.m.RHi  i  ) 

IF< INOFX.EO.llGO    TO    999 
NMAX=l 

""    fJfflf^WVnTIMMMm     999.999.79 
79     K'MAX  =  NMAX+1 

go  m  R9 

q90  CONTINUE 

|!B!f  r^/fvSOTff  ?J  IO.OlMI  +  1  I  1  .975.FN2.1  .05  .  / 
*(iMs'.0*(V'Er,Z**O.R)  ) )) 
,i4ll7fit3i.R43<.*V6G2*tni?**2)*U.+0.l*XI 
r,4=04*(  l  .-n.39fl»x ) 
n?=ui?*( 1.0/G*-l.o/G2l 
^2vir;irUl.HTl»n)-((ll.*mni«2)»Zl»RNl*Wl)/l«.»-(lll.*FTK1»».Ml* 

SUNlU',^rMTr4n?r;im!,.l(nn2U?),Z2*RNieW2»/l4..-MU.*FTK2*RN2. 

I    w?»( inRO?**?) -1002**2) I ) /l*.l ) 
0V1«*.»FNV1 /FRS1 

4«^M^HT'-RNl*nnl-RN1*l=Nl*FTKl*mF0J-nDlJi 
*S2-Z2*(HT2-RN2*QD2-RN2*FN2*FTK2*IOF02-OOZ)] 

GSl«6l/*Sl 

r,<;?=02/4<;2 

RFi  =  nvi*r,!;i/viM 

Fi-0    00?S2RR-(4.11S/    10.**9)  >*RF1 
F21;":o02';2^-!4al5/(in.**9))*RE2 

I.P1=SL1*W1 

lVFr,l='f,D*?:/n'i?.,533R*(0I1**2)) 

rn-T/^,"7^p/i1l**nU"^'0.01«.  +  .<1.975*FN,.1.05./ 

*h66.'o*iv=r,i**o.R) ) ) ) 

03.!7**3l.  R434*VEr,i*(0Il**2)*fl.+0.1*X) 

G3=03*(1.-0.39R*X> 

Ol=HAl*(1.0/r,3-l.0/Gl) 

RMi  =  r,3/r,i  „   ,. 

!F( IN0FX.EO.2)    GO    TO    *5 
IF(WTl.ME.WTU)     WT1=WT10LO 

,7    WtUwFT  5-<UFTS-WTl  )*l  (  1 .n-OFXP .01  )  )  / <RM1-DEXPC01) )  ) 
27    WT%^T"-(WT^-WT2)*(li;-nFXP<02))/<R^-nEXP(02))) 

IF(04RMWET5-KT5)-[1.0n-9).LF..0>    GO    TO    222 

WET5=WT5 

GO    TO    27 
222    WT3=WT1+(WT5-WT6)*RM1 

WT4=UT?-(WT5-WTft)*RM2 

IFIWT6-32.0)    767,221.221 
767    WT7=32.0 

WT55=WT7-(KT7-WT2)*((l.-nEXP(D2)l/IRM2-DEXPID2ll) 


80 


757  VFG1  =  MWT55-WT7)*VEG2*  1012**2)  )/(<WT55-WT6)*(DI  1**2)1 
^ri1:^?f"7^VAll:.S!5?r?*<0.01«>*<<1.975*FNl+1.05,/ 

*(  1«^.0=(VEG1**0.B)  I  )  ) 

03ii7)s*l\.«ft3*«ver,i*i  nn**2)*  <  i.+o.i*x) 

R3=O3*(l.-0.39R*X) 

Dl  =114  1*1  1.0/G3-1.0/G1) 

wf  l=WT5S-(WT55-UTl  ).(ll  .O-OFXPIOl  )  )/  CRM1-DEXP (01 )  )  ) 

Wt^iT»fc7i^^?i:-S!^ii/i.M&i»»io*iii 

iF(nABS(wr55-wT5)-(i.on-9). 1.6.0)  go  to  223 

WT55=WT5 

GO    TO    757 
223  WT3«WT1  +  (WT5-WTM*RM1 
"3  5t*.WI2-(WT5-WT7)*RM2 

221     lF(WTl.NE.WTll)     GO    TO    47 
05=04-03 
WVE01=VEG1 
WVEG2=VEG2 
WWT1=WT1 
WWT2=WT2 
WUT3=WT3 
WWT4=WT4 
WWT5=WT5 
WWTft=WTfc 
WWT7=WT7 
GPM11=GPMI 
II11=H1 
IIMI'MM 

RM11=R"1 
011=01 
011=01 
Gll=r,i 
022=02 
f,22=G2 
033=03 
G33=G3 
044=04 
G44=G4 
055=05 
GPH22=GPM2 
1122=1)2 
11622=1142 
PM22=R«2 

WTil=(0.5=[WWT1*WWT3)+4ftO.) 

WT42=<0.5*(  WWT2  +  WWT4)+4f>0.) 

45    CONTINUE  .       .....  „ 

IF(STl.NE.STU)     ST1  =  ST101_D 

30    STA-SET5-'("t5-St1  )*l  (  1 .0-OFXP  (01  I  )  /  IRM1-0EXP  (Dl  )  )  ) 
lT5«STfi-(STA-ST2)*( ( 1 .n-OEXP < 0? > ) / I RM2-0EXP < D2 ) ) ) 
!F(OARMSET5-ST5)-(1.00-9).LE.O)    GO    TO    225 
SET5=ST5 
GO    TO      30 
225    ST3  =  STl  +  (<>T5-STf.)*i«Kl 
ST4=ST2-(ST5-STfc)*RM2 
IF(STl.NE.STll)    GO    TO    47 
SVEG1=VEG1 


SVFG2=VEG2 

SST1=ST1 

SST2=ST2 

sst3=st3 

SST4  =  <;T4 
<,<;T5  =  ST5 
SST*,  =  ST6 

stai  =  io.5*(<;sti  +  <;<;t3i+<.60.> 

STA?=(0.5*(SST2+SST4)+AhO.) 
47    CONTINUE 
*  TAIL    HOURS    (!NOi=X.WTl,WT3.STl.ST?.N,OLnT.ni.nH.RNEWT.RNEWH,TAHGR. 

1P(  INDEX. 5.0. 2)  GO  TO  111 

1FIUT1. I.E. 70.0.  AND. WT1. I.E. WT2  1  GO  TO  999 

111  "lP(VTl.GE.50.0.AND.STl.GE.ST2l  r,o  TO  45 

113  continue 

nsft  =  n<;ft+i 

IF(NSCT-?I  997.997,1000 
1000  CONTINUE 

wofni=3o,sR5/WTai 

WnFN?  =  39.ftf>5/WTA2 

WS1  =  W0FN1 /">?.* 

"APo"=(Fi>*^pi*('-,<;l*«?)^(!OVl/SHl)»*0.'H«((S1.1/SHl  1**0. 61*27. ft923) 

^/pn22rF2irP^(^i^t:MnV?,/£H2)»0.4..<(S..2/SH2,»0.6.*27.6923) 
1     /(r)V?*wS/!*(,5,2?«i!  10.*"10M  ) 

wr,RO=   TAM"R*ooT/ionnriO. 
wr.Rr,.   TAHr.Hiiir.T/ionono. 

wnFci  =  (0H74/Unp>cu  srFw^wAPOi  )*wffi_/  f6356,0*P61 ) 
5nFC3=  0.7A6*0PKW  *r,FliF:sWAPn?)«WEFI./(6356.0*FE2) 
HRRl  =  1.0R*r.F«S»(WWT3-WWTl  1 

a^V-1l".M«"£*WvUl/»!o?«l.T5l*Ct.*0.1«l«IO.OOS»«l«ln». 

»  n'(l»7Ssigi-l  6  1+0. *)+  (  (  35.7R/6*.*)*!  WVEG1**2)  ) 

*  O.0R75*(W?-l. 01+0. 31+1  (35.7R/64.4)*(WVEC,2**2)  I 
WPFG-?  0*WPFG2 

WO  =  P=l0.7  +  Asr,PM2*WPEG*OPKW  1*WFFI./(3960*PE) 

SnFNl=39.ftfl5/STAl 

SnFN?  =  39.<,R5/STA2 

SS^SOFNl/f,?.* 

SAPni.lFli? Pi«(fi<;i**2)* I  (0V1/SH1 1**0.4 )*I(SL1/SH1I**0.6)*2T.6923) 
lSip^FHrP5?*"*S?**2)*Mn{r2/SH2>..0.*).USL2/SH2)..0.6).27.6«3l 
1-/(DV?*SS2*I5.2?*(10.**101 )) 

Sf,R  =    TACGP.snpTR/12000. 
SG°l=l  .0«sC  =  "S*l  SPT1-SST3) 

^OFrri(0.74+,*r)p"KW  *CFMC,-SAP01  )eSEFI./<*35f>.0*Ffl  I 

S0FC?=(0.746*nPKW  *r.FMFsSAPn?l*SFFl/  1 6 35A  .  0*FF2  I 

■  SPFGl=  l.*9S)*((  <!VFr,l/3. 01**1.75  )*  I  1.+0.1*X)*IO. 0039*«1*Z1*12+ 

*  0  0R75*(W1-  1.0  1+0.31  +  1  (3«.7R/64.A)*IUVFG1**2)  I 
%p4o?.(i.l^5)*|(SVFfi?/3.0)**1.751*(l.+0.1*X)*<0.0039*W2*22n2+ 

*  0.0875*  (W2-1.01 +0.3  1  +  1135. 7fl/64.M*(WVEG2**2)  I 


o2 


|KK1o!%UgI8!wM«WIW    "SEFL/.3960.PE. 

tgro=wrro+sgr 
tgrg=ugrg+sgr 
lffi!wnlc!^nEci*wneci+so6C2*M06P*soep 

TA0F=T0F*34  13  .0/0.025 

fnro=trro-toe 

fnrg=tgrg-toe  .,,,, 

FNRF=TGRF-TOE  „m*7  1  I  /  (  ( Wl«*0.333  )  »  "  'fl  /1000  .  )**0  .13  ' 1 

TrRR=nr.HT  *wr,Ri/looooo.O 
TC«C*r>r.CT*SRRJ/12000.0 

S&fiai^clirticweicPfBiCLrtiweiw 

A  =  A'l+A2  .  ,„„ 
TAGR=TAHOR+TACGR 
TANRR=TAGR-TAOE 

KW&lff£tof«i 

PRINT  1" 
,01  FtWW'l.'.20X,'  (WINTER  RESULT  )  •) 

t  WfiR2^GRn"r;RG.NGRE,W0eCl.WOF.C2.HWf7 

1FI05.GT.O.O)     GO    TO    107  .    . 

107    WRITEIA.10R) 
PRINT    1<5 

103   F0RM4T("*'.20X,'     (SUMMER    RESULT    )     •) 

PRINT    1«  „_„    rrT.     C<-Tc;    tu*    RPMl  .GPM2.01  ,02  .03.04. 

WRITE!  A.  10*1    |S":SST4.SST5.SSTA,GP«1    r,Pu|ji|il.jfi  f 

PRINT    1« 
PRINT     105 
105    FORMAT! ■»■ .20X,        

%  TGRPA.PAYRO.PAYRG.PAYRE.TANGR.A 

PRINT    1» 
1    F0RMATII3)       ^_ 
3    FORMAT (F7.1  ,FR.l) 

II  SBStfi::.':  «"S  •  7  &»iF»&™  "•'-. 
S^3.*7!.,i.,:8^35{IJSj|585 

34  FORMAT!  '*,.10X,  THIS  CASE  FOR 

cn  FORMAT  I  '*■  . ]0X, ■ TH1 S  CASE  hUK 

\1    FORMAT  •»'  10X. -THIS  CASE  FOR 

4Q  FORMAT!'*' .10X, 'THIS  CASF  FOR 

52  FORMAT! •*'.10X, 'THIS  CASE  FOR 


■(COMMON  AND  CONCLUSION  RESULTS)') 


AIR  VELOCITY  IV)=100  FPM  ' 

AIR  VELOCITY  IV) =200  FPM  ' 

AIR  VFLOC ITY  (V)=300  FPM  ' 

AIR  VELOCITY  IV)=400  FPM  ' 

AIR  VELOCITY  (V)=500  FPM  • 

AIR  VELOCITY  IVI-fcQO  FPM  I 

AIR  VELOCITY  (V)=700  FPM  ■ 


as 


54  FORMAT! '4> ,10X, 'THIS  CASE  FOR   AIR  VELOCITY  IV)=SOO  FPU  <) 


200  FORMAT! 

/.ox. ■ 

C^MS 

=  . 

.F12.5.5X,' 

CFME 

«■  . 

E12.5.5X.' 

WT1 

=  ' 

.E12.5, 

* 

/.ox. 

WT2 

.F12.S.SX, 

ST1 

F12.5.5X, 

?T2 

.E12.5, 

/  , OX  ,  ' 

V  -  -, 

.F12.5.5X,' 

X 

F17.5.5X, ' 

001 

•E12.5. 

* 

/.ox. 

on? 

s 

.F12.5.5X, 

VM 

=  i 

F12.5.5X, 

VA? 

= 

•F12.5. 

• 

/  .ox. 

FN  1 

= 

.E12.5.5X.' 

CIO? 

B  1  . 

E12 .5.5X, ' 

SHI 

=  ' 

•F12.5, 

4 

/•nx, 

SH? 

• 

.F12.5.5X, 

Sl.l 

s  * 

E12.5.5X, 

SL2 

= 

,F12.5, 

4 

/  .ox , 

RN  ) 

.F12.5.5X, ' 

R'J? 

E12.5.5X, • 

wj 

•F12.5, 

* 

/,QX, 

H? 

.F12.5-.5X, 

FEl 

E12.5.5X, 

FE2 

•E12.5, 

* 

/  .ox , 

PE 

- 

.E12.5.5X,' 

WEFI 

■  > 

E12.5.5X,  • 

SEFI. 

=  * 

,E12,5, 

* 

/.ox. 

F  ICP 

= 

.F12.5.5X, 

E  If.P 

-' 

F12.5.5X. 

EICL 

- 

,E12.5, 

4 

/  ,0X  , 

c  I  cg 

.F17..5.5X. 

ETCF 

E12.5.5X, 

TK1 

,F12.5, 

* 

/.ox, 

TK2 

= 

,  F  I  ?  .  5  ,  5  X  , 

OOT 

=  * 

E12.5.5X, 

OGT 

= 

,F12.5, 

4 

/  .ox  , 

OFT 

- 

.E17.5.5X,' 

opkw 

s  »  , 

E12.5.5X, 

OPTR 

= 

.F12.5, 

* 

/.ox. 

OCHT 

- 

.E12.5.5X, 

nr.r.T 

B  1 

E12.5.5X, 

VI  SI 

a 

,E12,5, 

4 

/  ,ox. 

VI  S2 

- 

.F12.5.5X. 

FTK1 

=  1 

E12.5.5X, 

FTK2 

= 

,E12.5) 

102  FOR 

4 

MAT! 

/,ox, 

T3 

~ 

.F12.5.5X, 

T4 

=  • 

F12.5.5X, 

T5 

= 

,E12.5, 

/  .ox  , 

T*, 

.F] 2.5.5X. 

r,PMi 

E12.5.5X, 

GPM2 

.E12.5, 

4 

/.ox. 

01 

- 

.F12.5.5X, 

02 

=  1 

E12.5.5X, 

03 

= 

,E12.5, 

4 

/  ,ox , 

04 

■E12.5.5X, 

05 

F12.5.5X. 

Gl 

.F12.5, 

a 

/.ox. 

<",2 

- 

,E12.5 ,5X, 

ft  3 

a  1 

E12.5.5X, 

r.4 

= 

,E12.5, 

* 

/  .ox. 

ni 

- 

,F1?.5 ,5X, 

02 

=  ' 

E12.5.5X, 

til 

= 

.E12.5. 

* 

/.ox. 

ii? 

,E12.5,5X. 

IIA1 

E  1  2  .  5  ,  5  X  , 

n»2 

,E12.5, 

4 

/  .ox , 

RMl 

.F12.5.5X, 

RM? 

E12.5.5X, 

VEG1 

.F12.5. 

* 

/.ox, 

VEG? 

= 

•F12.5.5X, 

PEfil 

e  » 

E12. 5,5X , 

PEG? 

= 

,E1?.5, 

* 

/  ,ox , 

DFr, 

- 

.E12.5.5X, 

AP01 

s:  ' 

E12.5.5X, 

AP02 

= 

.E12.5, 

4 

/,  ox, 

or-p 

,FI 2 .5 ,5X, 

RSI 

E12.5,5X, 

RR2 

,F12.5, 

* 

/  .ox , 

ORO 

.E12.5.5X, 

ORG 

E12.5.5X, 

GRE 

,E12.5, 

* 

/.ox, 

OFC  1 

= 

.F12.5.5X, 

0FC2 

s  ' 

E12.5.5X, 

T7 

= 

■F12.5) 

10B  FORMAT  1 

* '  »BX  , 

(THE 

Rl 

-PASS  VAI.VS 

IS  01^ 

1  IN  WINTE! 

OPERA' 

ION) ■ ) 

10*  FORMAT! 

/.ox. 

T3 

- 

,F 1?  .5  ,5X, 

T4 

=  ' 

,E12.5,5X, 

TS 

= 

,E12.5, 

4 

/  ,QX  , 

Jh 

.E12.5.5X, 

GPM1 

E12.5.5X, 

GPM2 

.F12.5, 

/  ,  0  X  , 

01 

- 

.E12.S ,5X, 

02 

=  ' 

,E12.5,SX, 

03 

= 

'  .E12.5, 

4 

/  .ox  , 

04 

•  ■=)?. 5, 5X, 

0FN1 

F12.5.5X, 

Gl 

.E12.5, 

* 

/.ox, 

G2 

- 

.F12.5.5X, 

G3 

s  * 

.E12.5.5X, 

'G4 

a 

'  ,F12.5, 

4 

/  .ox , 

D! 

,E1?.5,5X, 

02 

E12.5.5X, 

Ul 

,E12.5, 

4 

/.ox, 

II? 

B 

.F12.5.5X, 

DAI 

=■  1 

.  F12.5.5X, 

•IJA? 

= 

'  ,E12.5, 

4 

/  .ox. 

RMl 

- 

.F.12.5.5X. 

RM2 

s  ' 

E12.5.5X, 

VEG1 

-- 

,F12.5, 

* 

/,ox, 

VEG2 

- 

,F1? .5,5X, 

PEG1 

B  * 

■E12.5.5X, 

'PEG? 

= 

'  ,F12.5, 

4 

/  ,ox  , 

per. 

.E12.5.5X, 

AP01 

E12.5.5X, 

AP02 

,E12.5, 

4 

/,ox. 

'OFP 

- 

.F12.5.5X, 

GR1 

=  1 

.E12.5.5X, 

'GR2 

= 

'  ,E12.5, 

4 

/.°x. 

OECl 

- 

.F12.5.5X, 

0FC2 

=  • 

E12.5.5X, 

GR 

- 

.F12.5) 

106  FORMAT! 

/.ox. 

•HT1 

s 

.F12.5.5X, 

•HT? 

--  ' 

•E12.5.5X, 

'  HOI 

- 

1  ,E12.5, 

4 

/.ox, 

WO? 

,612, 5, 5X, 

21 

El? .5  ,5X, 

17. 

,F12 .5, 

/.ox. 

'TGRO 

•  ,E1 2.5.5X, 

'TGRG 

,E12.5,5X, 

'TORE 

' ,E12.5, 

4 

/,ox  , 

TAOE 

= 

.E12.5.5X, 

FMRO 

-    ' 

E12.5.5X, 

ENRG 

- 

■E12.5, 

4 

/.ox. 

>FNRF 

- 

■  ,F1?.5,5X, 

•  Tr.RB 

=  » 

,E12.5,5X, 

'TCRC 

- 

' .F12.5, 

4 

/  ,ox  , 

TCER 

= 

,F1?.5,5X, 

FIT. 

=  ' 

E12.5.5X, 

TAHGR= 

.F12.5, 

4 

/.ox. 

'TACT, 

R  = 

• ,F1?.5,5X, 

•  TAGR 

a  ' 

.E12.5.5X, 

' TGRPA= 

' ,F12.5, 

* 

/  .ox  , 

PAYRO= 

.E12.5.5X, 

PAYBG=' 

E12.5.5X, 

PAYBE= 

,E12.5, 

4 

/,ox, 

■TANGR= 

1 .E12.5.5X, 

•A 

=  ' 

•E12.5) 

MM  = 

MM  +  1 

IF 

MM-KK)  5,5, 

7 

7  CONTINUE 

STOP 

C      i«s*#satt**:*a* 

4te£* 

44 

fcfi444444444 

»«tt 

^  =  i 

4444444444 

44444 

*ft 

44444444 

C      4 

THIS 

S'IRROII 

r  IN; 

HO 

IRS  IISEO  FO 

(  TABU'.A 

IING  NEW  A 

?RAY 

3F 

TEMPS. 

C         4 

AND 

ARRAY  0 

F  COR 

RE 

SP0ND1NG  HO 

JRS. 
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&#*##*=*#=**#**=**** 


e***#*ft*#**#**#************************** 


SURROUTINE    HOURS! INDEX, WTl.WT3,STl,ST3,N,OLOT,OLOH,RNeWT,RN6WH, 
^SURRUUIINS    MulK^llT4HGi,iTicr^,CFMS,,_L,i.|_u,wT  101.0.  ST10LD.NM1X) 

OI,«FNSinNRn|4riT?5o7r6LnH(50),RNEWT(50),RNEWH(50) 

IF( JNnFX.E0.2l  go  TO  31 

THIN=WT1 

TM1VX  =  WT3 

GO    TO    32 

31  TMIN=ST3 
TM»X=ST1 

32  CONTINUE 
N<=N-l 
00  200  1=1. NX 

"ilF1'(TMIN.(-,F.OinT(T)  .SNn. THIN. I.E. 0I.0T(!+1  )     )    60    TO    14 
IF     (TMtX.r.FlniOTIJ-n.lND.THax.I.E.ni.DTIJ))     GO    TO    15 
GO    TO    200 
l*    hm'nL(OI  DHtKKl  +  1  )-ni  r>H<KKl)  )/ IDIOT (KK1+1)-0L0T(KK1)  )* 
I     rTM]M— HI  0TIKX1)  l+OI.DHIKKl) 
"     J-l 

1  (TMJx-OLnt (Kk2) i+nCnH(KK2l 

200  CONTINUE 

00    100    1=1. NX 

J=I*1 

|F[M'!!lI>N!rr:??ni?OH??)Is,"o!HMiM.i.E.Oi.nH(l  +  l)l    GO   TO   40 

IFIHMix.r.F.OLDHI  J-l  )  .4Nn.HMax.LF.0L0H{J  )  )    GO    TO    30 

GO    TO    100 
41    ^FIhJiNJ  e.ni.0Hm.4NO.HMlN.R6.Ol.nH(I+l))    GO    TO    40 

IFIHM1X    CF;oCDH(J-l).iNO.HM»X.GF.OLOH(JM    GO    TO    30 

GO    TO    100 
AO    K1=I 
30    <?=J-l 
100    CONTINUE 

NN  =  i<2-Kl 

NN]=NN+1 

NN2=NN+2 

RNFWHI 1  I =  HMJN 

RN=WH(NN2)  =  HMAX 

RNFWTI1 l=THTN 

RNFWT (NN2) =TMAX 

00     RO     L=2.NN1 

RMFWH(I    )=Tl.nH(Kl  +  l.-l  ) 

RNFWTd   )=0L0T(K1+L-1) 
RO    CONTINUE 

l"F(  INOFX.C0.2I    GO    Tn    95 
HSRi(l,0«*CFHS)*(THAX-RNEWTIU))*RN6MHIt) 

WT1=RNFWT1I.  +  1) 
WT10LO=WT1 
T6Hr,^  =  TlHGR*HGR 
U.=IL+1 
GO    10    17 

05    rGR=1.0R='FMS*(RNFWT(NN2-J)-THIN)*RNEWHINN2-J) 
STl=OLDT(NMtX-LLL) 


BE: 


STiOi.O'STl 

Tarr,R=T4Cr,R+CGR 
i  i  i  »UJ.+l 
Q7    CONTINUE 
RETURN 
END 
SENTRY 
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Appendix  C 

Determination  of  VS 

The  equations  of  the  curves  in  figure  28  were  determined  by  the 
computer  program  (shown  in  Appendix  D)  as  shown 

VS  ■  (298.9753  -  11.59322  x  10"3  H  +  0.1449816  x  10"6  H2)  for  AV  =  300 

VS  =  (399.9628  -  15.44916  x  10"3  H  +  0.2769957  x  10"6  H2)  for  AV  ■  400 

VS  =  (499.7739  -  18.673  x  10"3  H  +  0.192449  x  10"6  H2)  for  AV  =  500 

VS  =  (600.0673  -  22.18371  x  10"3  H  +  0.250031   x  10"6  H2)  for  AV  =  600 

VS  =  (699.8415  -  25.33304  x  10"3  H  +  0.2520391  x  10"6  H2)  for  AV  =  700 

VS  ■  (800.0952  -  29.20782  x  TO-3  H  +  0.336966  x  10"6  H2)  for  AV  =  800 

VS  =  (896.5751  -  27.84579  x  10"3  H  -  0.219282  x  10"6  H2)  for  AV  =  900 

But  since  VS  =  f(H,AV),  the  above  coefficients  were  used  as  input 
in  the  same  computer  program,   and  equation  30  was  found.     With  H  and  AV 
(input  data)  known,  VS  is  determined. 

A  smaller  value  of  VS  corresponds  to  smaller  air  flow  rate  and  a  large 
value  of  VS  corresponds  to  a  larger  air  flow  rate.     VS  was  calculated  by 
equation  30. 

Determination  of  R 
The  equations  of  the  curves  in  figure  27  were  determined  by  a  computer 
program  (shown  in  Appendix  D)  as  shown. 

R  =  (0.75771719  -  0.33494178  x  10"3  VS  +  0.13688623  x  10"5  VS2)     for  G=1.0 
R  =   (0.77299494  -  0.35435311  x  10"3  VS  +  0.16081322  x  10"6  VS2)     for  S-1.1 


S7 


BS 


R  =  (0.79764146  -  0.3404084  x  10"3  VS  +  0.12507371  x  10"6  VS2)  for  G=1.2 

R  ■  (0.81985146  -  0.37411135     10"3  VS  +  0.14926809  x  10"6  VS2)  for  G=1.3 

R  =  (0.83567923  -  0.38196915  x  10"3  VS  +  0.15071691  x  10"6  VS2)  for  G=1.4 

R  =  (0.84511191   -  0.38216077  x  10"3  VS  +  0.14893794  x  10'6  VS2)  for  G=1.5 

R  =  (0.86138093  -  0.40705432  x  10"3  VS  +  0.16623653  x  10"6  VS2)  for  G-1.6 

R  =  (0.88449264  -  0.44027995  x  10"3  VS  +  0.18759965  x  10"6  VS2)  for  G=l .8 

R  =  (0.91029024  -  0.48727333  x  10"3  VS  +  0.220064  x  10"6  VS2)  for  G=2 

Since  R  =  f(VS,  G),  the  coefficients  in  the  above  equations  were  used  as 
input  in  the  computer  program  (see  Appendix  D),  and  equation  31  was  found. 
Parameter  VS  was  taken  from  equation  30,  and  parameter  G  can  be  determined 

as 

P  _  Larger  air  flow    , 
Smaller  air  flow  - 

The  smaller  value  of  VS  with  G  value  were  substituted  in  equation  31,  and 

R  was  found, 
a 

Then, 

R.  =  R  /G 
b    a 

Determination  of  APD 

The  computer  program  shown  in  Appendix  D  was  used  to  find  the  equation 
of  APD  for  a   specific  value  of  H,  by  selecting  three  arbitrary  points  on 
each  curve  in  figure  29. 

APD  =  (-.7787674  x  10"1  +  0,830931  x  10"3  VS  +  0.1757019  x  10"5  VS2) 

for  H  =  0  ft 
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APD  =  (-.1104487  +  0.1141962  x  10"3  VS  +  0.1934748  x  10"5  VS  ) 


for  H  =  5  x  103  ft 


APD  ■  (-.8850657  x  10"1  +  0.1127086  x  10-3  VS  +  0.2726474  x  10"5  VS2) 


for  H  =  104  ft 


As 

APD  =  f(H,  VS),  the  coefficients  in  the  above  equation  were  used  as 
input  in  the  computer  program  (see  Appendix  D),  and  equation  32  was  found. 


Determination  of  F 


Table  1  represents  the  relationship  between  Fs  and  H.  Three  arbitrary 
values  were  selected  from  the  table  to  form  three  simultaneous  equations.  These 
equations  were  solved  for  the  coefficients.  The  equation  of  F$  becomes  as 
in  (33). 


APPENDIX  D 
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c  «»♦**,*»**•«««***»»••«    pv    iFa(.T    sf)uiRE    P0LYN0MIiL    METHOD)-  « 

r  «  N=     THF     NUMBER     OF     OATA     POINTS.  CTTT  =  n  • 

{  I      7l=  *TVMNnoro?p"nW§.Dftt*o?SKT«lirrtSf«KS  to  be  f.tteo.  j 

c  •  "ttl:  Yhf  BSVgWtS  \--th  8«4  p°o'inTt:   ,f   all  po.nts  are  . 

C  KmENSION   M3O.301.BI  30.30). HI30 )  ,4C(30,*J,m3O)  ,CC30, 301 

MM  =  0 

RFsrns.ioo)  li.  __  .„„„ 
5  lF|MM.r,F.i  I  )  <-,n  td  loon 
10  RFin   i 5. ?on>   n.m.l 

IF(MM.mF.O)    RO    TO    25 
15    WRITE! ft. '00) 

RMrVdTinni      H(i), win.!    ACU,J),J=1.U 

20  WRITE  '(f.U  00 1   H(!),wm.  ACH.J) 
f         FOLLOW  11  STATFMENTS  TO  FORM  THE  COEFFICIENT  MATRIX  OF  THE 
C      '.'.    NnnMAl.  EQUATIONS 
C 

25  00  30  1=1. N 
30  CI  1,1 1=1.0 
MP1=M+1 
00  35  .1  =  2, MP1 
00  35  1=1. N 
35  C(I,.I)=C(I.J-D*  HI  I) 

oo  4n  i=i ,mpi 
no  40  j=i.mpi 
ai  i  ,.i)=n.o 

40  Mltj)=M1I?J)+C(<.!)*C(K.J)*W(K) 
P  FOMOWINr,  5  STATEMENTS  TO  FORM  THE  CONSTANT  MATRIX  OF  THE 

C       '.'.     NORMAL  EQUATIONS 
C 

00  45  .1=1, L 

no  45  i=l ,mpi 

R(  1  ..11=0.0 

45  Rutj)=B,(l'i!j)+C(K,I)*AC(K.J)*W(K) 

C      ..  CALL  SUBROUTINE  C"APR  TO  SOLVE  THE  NORMAL  EQUATIONS 

CALL  CHAPR(A,MP1,B,L.0ET) 
WRITE  (ft  ,ft00)N 
WRITF(ft,700)  M 
00  50  J=l,L 
00  50  1=1, MP1 

50  WRTTE(f.,800)  II.B(I.J) 
WRITElft.SOO) 
MM=MM+1 
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c 


GO  Tn  5 
inn  FORMAT! !A) 

?nn  cnwuuiii 

300  format(?fio.a,fio.5!  _     ..,.,,. 

4nn  cnRMiTiM'.i'x.'H'.pnx.'WEir.HT'ty^x.'AC1//) 

500    FnRMAT(RX,F12.*.12X.Fl?.ft.nx,|15.6/> 

hnn   format    (///.'    niimrer   of   given  nATA   points   «    '.12) 
7nn   format(/7x. 'Ofgref   of   polynomial   =    ',12.///),, 
Ron  pormat i5x,i2.>    decree   coefficient  =    ".eir.r./i 
qoo   FORMAT!/, ioh   -   -   -   -   -,///) 
1000    STOP 
END 


s^nonN|  cmapp  ia.mpi.^.l.oet^^^^  jNn  STMJU#t|Me„  Ens. 

niMFNMn"  HI30,m'l  ,R< 30,30)  , I PVnT! 30  I, INDEX (30, 2), PIVOT  130) 

COMMON  1PV0T , INDEX , PIVOT 

EDO  I  VALENCE  ( IROW, JROUl , ( ICOL.JCOL) 

C      FOLLOWING  3  STATEMENTS  FOR  INITIALIZATION 


C 


57  DET=1. 

DO  17  .1  =  1, MP1 
17  lPVnT(.I)=0 

00  135  1=1, MPl 


C      FOIIOWING  12  STATEMENTS  FOR  SEARCH  FOR  PIVOT  ELEMENT 

c    t=0; 

no  o  .1=1  ,mpi 

IF  I  JPVOTI  .M  .EO.l )  GO  TO  9 
13  00  ?3  K=1,MP1 

IF  (IPVOT!KI-l)  43.23.R1 
43  IFIARS1T I .GE.ARS! A! J,KI ) >  GO  TO  23 
R3  IROW=.! 

!COL=K 

T=A(J,K1 

23  CONTIN'IE 

9  CONTINUE 

IPVOT! ICOL)=IPVOT( IC0L1+1 

C      FOLLOWING  15  STATEMENTS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

c 

IF  | IROW.EO. ICOL)  GO  TO  109 
73  DET=-Or:T 

DO  I?  M=1,MP1 

T=A( IROW.M ) 

A( IROW,M)=A( ICOL.M) 
12  A(ICOl..M)=T 

IF  (L.I.E.O)  GO  TO  109 
33  no  2  M  =  l  ,1 

T=R( IROW.M) 

RtlROW  ,M)=R(ICOL,M) 
2  R(  ICOI.,M)=T 
100  INDEX!  1. 1) =IROW 

INnEX! I ,2)=IC0L 

PIVOT ( 1)=A( ICOL, ICOL) 

DET=DET*PIVOT( I ) 

C      FOLLOWING  6  STATEMENTS  TO  DIVIDE  PIVOT  ROW  RY  PIVOT  ELEMENT 
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M  I  cm..  icoi.)  =  l. 
nn  ?05  m=i ,mp  i 

205  M!CnL.«)=MICOI.,M)/PIVnT(I) 
IF(1..l.F.0>  GO  TO  347 

52  S?IC0l."«)=MIC0L.«)  /PIVOT!  I) 

C      FOLLOWING  10  STATEMENTS  TO  REDUCE  NON-PIVOT  ROWS 

°  347  00  135  1.1  =  1. HP1 

I F( 1 I.EO.ICOL)  CD  TO  135 
21  T«M1.f.lC0U 

AIL1 . ICOL)=0. 

00  R«  M=l.WPl         

RQ  »IL!.Hl«»(Lt.M)-»( ICOU.H)*T 

IP(I  .1  F.O)  CO  TO  135 

in  no  as  m=i,l      „,„,_ 

6«  R(l  I  ,M)=«(L1  ,M)-B(ICOI..M)*T 
135  CONTINUE 
C      FOLLOWING  11  STATEMENTS  TO  INTERCHANGE  COLUMNS 
222  00  3  1=1. MP1 

IFMMNnEX(K,l).E0.IN0EX(M.2)i  GO  TO  3 
1Q  .|ROW=INnFK|M,l  | 
JC01  =IN0EX(M,2) 
00  54=)  K=l  ,MP1 
T  =  A(K,.IROW> 
A(K,.l*OW)=A(K,.ICOL> 
AIK.JCOI.  )=T 
540  CONTINUE 
3  CONTINUE 
Rl  RETURN 
ENO 
SENTRY 


APPENDIX  E 


94 


95 


r  *##****♦***#»***♦************«************************»*********** 

C  »              -((   O-PIPE  THERMAL  RECOVRY  UNIT  ))-  * 

C  *     H=  JOR  SITE  ELEVATION  IN  FEET  * 

C  *  CFMS=  SUPPI  Y  AIR  FLOW  IN  CFM  * 

C  *  CE"E=  EXHAUST  AIR  EI.OW    IN  CFM  * 

C  *    FA=  FACE  ARFA  (ONE  SIDE  DF  PARTITION)  IN  SO. FT.  * 

C  »   AVS  =  ACTIIAI   AIR  SUPPLY  VELOCITY  IN  FPM.  » 

C  *   AVE*  ACTIIAI  AIR  EXHAUST  VFI.OCITY  IN  FPM.  * 

C  *     G*  MASS  Fl  OW  RATIO  OF  THE  TWO  AIR  STREAMS  (LARGER/SMALLER)   * 

r  *   VS1*  STANDARD  AIR  VFIOCITY  nF  SMALLER  CFM  * 

C  *   VS?*  STAMnARD  AIR  VFIOCITY  OF  LARGER  CFM  * 

C  *    RA*  RFCDVFRY  FACTOR  OF  SMALLER  AIR  FLOW  » 

C  *    RR*  RECOVERY  FACTOR  OF  I  ARr.FR  AIR  FLOW  * 

r  *  APD1=  ACTUAL  PRESSURE  DROR  FOR  SMALLER  AIR  FLOW   IN  INCH  W.G.   * 

C  »  AP02=  ACTIIAI   PRESSURE  nROO  FOR    I  ARGER  AIR  FLOW    IN  INCH  W.G.   » 

C  «    FS=  AIR  FLOW  CORRECTION  FACTORARGER  AIR  FLOW  IN  INCH  W.G.     » 

r  »    wf  =  WJNTFR  TFMPFRATIIRF  IN  DEGREE  FAH.  * 

C  *    ST  =  SUMMFR  TFMPFRATIIRE  IN  DEGREE  FAH.  * 

C  *    FF=  FAN  EFFICIENCY  IN  (2)  * 

C  *  FICll"  INITIAL  COST  OF  THE  O-PIPE  THERMAL  RECOVERY  UNIT  IN  IS)   * 

C  *  F'ICI  ■  I  AROR  INITIAL  COST  IN  (I)  * 

C  *  EICF=  FANS  INITIAL  COST  IN  IS)  * 

C  •    GR  =  ENERGY  GROSS  RECDVERD  IN  RTIIH.  * 

C  «CFMSI=  AIR  SUPPI  Y  BY-PASSED  IN  CFM,  * 

C  *  TGRO*  TOTAI  GROSS  RFCOVFRED  1011.1  IN  IS/YR.)  * 

C  *  TGRG=  TOTAL  GROSS  RFr.nvFRFD  (GAS)  IN  IS/YR.)  * 

C  *  TGRF=  TOTAI  GROSS  RECOVERED  (EI.EC.)  IN  IS/YR.)  * 

C  *   TOF*  TOTAI  OPERATING  FNERGY  COST  IN  (S/YR.)  * 

C  *  TAOF*  TOTAI   ANNUA'  OPERATING  ENERGY  IN  IRTU/YR.)  * 

C  *  FNRO=  NFT  RECOVERY  (Oil.)   IN  (S/YR.)  » 

C  *    FNRG*  NET  RECOVERY  (GAS)  IN  (S/YR.)  * 

C  *  FNRE*  NET  RECOVERY  (ELFC.)  IN  IS/YR.)  * 

C  *  TCRR*  ROIIFR  CAPITAl   RFDIICION  IN  It)  * 

C  *  TCRC»  CHII  I  ER  CAPITAl   REDHC10N  IN  (SI  * 

C  »  UFEl  =  WINTER  FOUIVAI  FMT  HEATING  LOAD  IN  IHRS./YR.)  * 

C  *  SFFL=  SUMMFR  FOUIVAI  ENT  COOIING  LOAD  IN  (HRS./YR.)  * 

C  *   DOT*  Oil  HEATING  COST  IN  IS/THERM)  * 

C  *   DGT=  GAS  HEATING  COST  IN  (S/THFRM)  * 

C  *   OET*  Fi  FC.  HEATING  COST  IN  (S/THERM)  * 

C  *  DPKW*  OPERATING  FMFRGY  COST   IFLEC.)   IN  (I/KW-HR.)  • 

C  *  DPTR  =  OPERATING  FMFRGY  COST  (El  EC.)  IN  IS/TON  OF  REFRIGERATION)* 

C  *  DCHT»  BOILER  CAPITAL  FOIIIPMFMT  RFDUCT10N  IN  I  S  /THERM  )  1GERAT  ION )» 

C  *  OCCT*  CHll'ER  CAPITAl  eoiiipmFNT  REDUCTION  IN  (S/THERM)  » 

C  *    0F=  OPFRATING  FNERGY  COST  IN  ($)  * 

C  *  TCER*  TOTAI  CAPITAL  EQUIPMENT  REDUCTION  IN  (S)  * 

C  »   EIC=  TOTAL  INTTIAI  COST  IN  Is)  * 

C  *TAHGR=  TOTAI   ANMHAI   HEATING  GROSS  RECOVERY  IN  (RTU/YR.)  » 

r  *TACGR=  TOTAI"  ANNUAI  COOIING  GROSS  RECOVERY  IN  (8TU/YR.)  * 

C  *  TAGR*  TOTAf  ANNIIAl"  GROSS  RECOVERY  IN  IRTU/YR.)  * 

C  sTANGR*  TOTAL  ANNUAL  NET  GROSS  RECOVERY  IN  (RTU/YR.)  * 

C  SPAYRO=  PAYBACK  (IN  CASF  OIL  USED  FOR  HEATING)   IN  IYRS.)  * 

C  *PAYRG=  PAYRACK  (IN  CASE  GAS  USED  FOR  HEATING)  IN  IYRS.)  * 

C  *PAYRE=  PAYBACK  (IN  CASE  ELEC.  USED  FOR  HEATING)  IN  (YRS.)  * 

C  *  *******  *******************************  **************  ************* 

IMPIICIT  R  =  AI  *R(  A-H.O-Z  ) 

DIMENSION  01  DT(50) ,OLDH(50)  .RNEWT ( 50 ) .RNEWH ( 50 ) 

REAL  WT11/0.0/.ST11/95.0/ 

MM*1 

READIS.A't)  KK 
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REACH  5, *>     H,CFMS,CFME.WT2.ST2.FEl.FE2,WEFl,SEFl,EICU, 

5  RFA0I5,*)  V 
WTI=WT11 

fFiMMlFn.i.AND.MM.i.T.2)  go  Tn  n 

IF(MM.FO.?.AND.MM.LT.3>  r,n  TO  9 

IF  (mm  .F0.3.AND.MM.I  T.4I  GO  TO  10 

1F(MM.F0.4.AND.MM.I.T.5I  GO  TO  11 

IFIMM.FO.5  .AND. MM. LT. 6  )  GO  TO  13 

IF(MM.Fn.6.ANf",.MM.1.T.7>  GO  Tn  46 

IF  (MM.F0.7.AND.MM.1.T.B  )  GO  TO  4R 

TF(MM.E0.R.AND.MM.LT.9)  GO  TO  53 
R  WRITF.(6,2B  ) 

GO  TO  35 
o  WR1TFI6.29) 

GO  TO  35 

10  WR1TEI6.34) 
GO  TO  3  5 

11  WRITE(6,50) 
GO  TO  35 

13  WRITF.16,51) 

GO  TO  35 
46  WRITEI6.49) 

GO  TO  3  5 
4R  WRITFI6.52) 

GO  TO  35 
53  WRTTF(6,54) 
35  PRINT  19 

PRINT  22 

W»1TF(6?2Q9)  H.CFMS.r.FMF.WTl.STl.WT2.ST2.FEl.FE2.WEFI  .SEFL.EICU, 
C"R"''6''   '  EiCLtlicF.DOT.OGT.DET.DPKW.DPTR.OCHT.OCCT.V 

CFMSI=0.0 

TAHGR=0.0 

TACOR-0.0 

NSET=1 

'.1-=1 

I.I  I  =1 
997  RFADI5.1)  1N0EX 

READ  (5,1 )  N 

00  12  I=1.N 
12  RFAD  (5.3)  OLDTI I) ,0LDH(I) 

CFHSU'CFMS 

IFICFMF.GT.CFMSll)    GO    TO    700 

F4=CFMSll/V 

G=CFMS11/CFME 

GX=G+l.0D-9 

M=GX 

GO    TO    701 

700  FA=CFME/V 
G=CFMF/CFMS11 
GX=G+1.0D-9 

701  AVS=CFMS11/FA 

tXm- (-7/ 30421  fc+1.0?7933»AVS-O.OOn0253223R*»VS**2)  + I  0.  0052771 19- 
i  o  SiZ-iZlii il  n*«  (  -4  I  SAV?*n.?R04574*l  0"  (  -7  )  *AVS**2 )  *H*  1-0.6582271* 
^6U"6K0.353UU*10«(-8)«AVS-0.3220021*10«*I-11)*AVS»*(2»)* 
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VS02  = (-7.30421 6+1.027933*AVE-0.O0OO253223R*AVE**2)-K 0.0052771 19- 
10  °327«h3*10**( -4) *AVE+0.2R04574»10*» 1-7)* AVE**2)»H* I -0.6582271* 
1?10*i(-(S)+0.3537444*10**(-R)*AVE-0.3220021*10**l-ll)*AVE**(2)(* 

3H**(2 ) 
IF  (AVS.GT.AVE)  GO  TO  702 
VS1=VS01 
VS2=VS02 
GO  TO  703 

702  VS1=VS02 

703  R"°o!2Il90998+0.336R77R*fi-p.pMR9942.(G**2))  +  «-0.366823/a 

r  t    l i^q?4?/l 1  a**"*)  )*G-(  .R7  691 / I 10**4    )*(G**2) )-vSl  +  < 1 < •£Hbuyc{/ 
r.(io.**A))-(  0.2522431/  (  10.**6  )  )*G+(  .  1094314/  I  10.  **6  1  >  *  <G**2  >  )  >* 
C(VS1**2) 
RR=RA/G 


IF( INDFX.EO.l)     GO    TO    999 

NKAX=l 
R9    CONTINUE  _       __ 

IFISTll-OI.OTINMAX))     999.999,79 
79    NM&X=NMAX+1 

GO    TO    R9 
999    |F(  INDC-X.E0.2)     GO    TO    45 

IF(WTI.N=.WT11)    WTl«WTinU> 

lFICFMSl  1.1  F.CFMFI  GO  TO  704 

WT3=WT1+(RR*(WT?-WT1)) 

WT4=WT?-(RA*(WT2-WT1> I 

GO  TO  7  05 

704  WT3=WT1 * (R6*(WT2-WT1) I 
WT4=UT2-(^f>*IWT2-WTl)  ) 

705  IFIWT4-32.0)  721,722,722 

721  WT4=3?.0 
CFMS=CFME*(WT2-WT4)/(WT3-WT1) 
rcMsi=cFM<;ii-CF^5 

WT5=(  (CFMS*WT3)  +  tf-FHSl*WTll  1/CFMS11 

722  IFIWT1.NE.WT11 )  GO  TO  47 
HWT1=WT1 

HWT2=WT2 
WWT3=WT3 
WHT4=WT4 
WWT5=WT5 
WRA=RA 
WRR=RR 
r.FMSS=r.FMS 
CFMS01=CFMS1 
AVS1=AVS 
45  IF(STl.NE.STll)  ST1=ST10LD 
lc(CFM<;.l  E  .CFHE)  GO  TO  706 

st3=st1+(rr*(st2-st1) ) 
st4=<;t?-ira*ist2-sti)) 

GO  TO  707 

706  ST3=ST1+(*A*(ST2-ST1)) 
ST4=ST2-(RB*(ST2-ST1) I 

707  CONTINUE 


IF(STl.NE.STll)    GO    TO    47 

SST1=ST1 

SST2=ST2 

SST3=ST3 

SST4=ST4 

SRA=RA 

SRR=RB 

&VS2=AVS 

cfmse=cfms 

JFF!>JTTLF.E7°O^.Sn:STr.iE.«T2,    GO    10-999 

1U    (P(sflMp.50.0.4ND.STl.G6.ST2)    GO   TO   999 

113    CONTINUE 

NSET=NSET+1 

1F(NSET-?I    99R,99RilOOO 
99B    CFMS=CFMS11 
GO    TO    997 
1000    CONTINUE 

WGPO=TaH,-,R*0OT/10000O. 
WGRG=TAHGR*OGT/ 10  0000. 
WGRE=TlHGRanET/ 100000. 
<;r.R  =  T4r.r.R*nPTR/12000. 

W0Fl'=m:74^n"KW*CFMJ*6"ni,)*WFFl./(ft356.*FEl) 
wr,f;!>-(0.746«nPKW«=CFME*SPn2)*WEFL/lh356.*FE2! 

714    wnFl1(oi74A«nPKW*CF>1<;*SP02)*WFFI./  Ift356. *FE1  ) 

unF?-iri ,  7'-f.*nDKw»r>MF*»pni)*wEFi./(A356.*FE2) 

710  WGPi  =  1.0R*r.FMSS*FS*(WWT3-WWTl  ) 
Wr,R?=l.0RarF«c*BS*(WWT?-WWT4» 
SGRl  =  1.0R#r.FMS5*FS*(SSTl-SST3] 
SGR2=l.0R*CFM=*FS*<SST4-SST2» 

S0E2=(0l746*DPKW*CFME*lPD2)«SEFL/l6356.*FE2) 

711  SOFTl  =  (0l.'4f,*OPKW*CFM<;S*4PO2)*SEFL/<f.35f,.eF|l) 
Ml    |nF2=(fl.74f,*nPKW*CFMF»4PDl)*SEFL/l6356.»F62) 

712  tGRii=WGRO+SGR 
TGRG  =  wr.Rr,  +  SGR 
TGRE«WGRE+SGR 
T0F=W0F1+S0E 1+WOE2+SOE2 
FNRO=TGR0-T0E 
FNRG=TGRG-TOE 

FNRE  =  Tr,RE-TOE 

TCRR=OCHT    twr.Rl/100000.0 

Tr.RC  =  OCCT»SGRl/  12000.0 

TCER=TCRB+TCRC 

EIC=EICU  +  EICI.  +  EICF 

TiOE=TOF*3413./0.025 

TAGR=T1HGR+T4CGR 

T»NGR»TAG*-T*nl= 

PAYP,n=:  (  E  I  C-TCER  > /FNRO 

PAYRG=(EIC-TCER)/ENRG 
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PAYBE=[EIC-TCER)/ENRE 
PRINT     10 

101    FORMAT  P*',20X,'     (WINTER    RESULT     )•) 

IJmii'ini]    UWT3    WWTi.UWT5.4VSl  .  AVE  ,  VS1  ,  VS?  ,WGR0  .WGRG  ,WGRE  .  APD 
1FICFMS1.0T.0.0)    GO    TO    723 

^U  fo«M»T?f*'.10X,     •     (THE    RY-PASS    IS    ON    )     •) 

PRINT    10 

PRINT    103 
103    FORMAT! •»»t20X,'        (SUMMER    RESULT)        •) 

WRITE(I^104)     SST3.SST4.AVS2.AVE.FA.VS1 . VS2 . SRR ,SGR 1 , SGR2, 
^WRITEIA.IO     I    ^Dj>;^n2;SRt,SRB,SnEl.SOE2.G.CFMSE.CEME 

PRINT    10 

105  FORMATS*', 20X,  •  (COMMON  AND  CONCLUSION  RESULTS)') 

WRIT  Ell" 106  I  TGRO.TGRG.TGRF.TOF.TAOE.ENRO.ENRG.ENRE.TCRB.TCRC. 
#WRITE(6.106I  ^'^^i'c^AHGR.TACGR.TAGR.TANGR.PAYRO.PAYRG.PAYBE 

PRINT  1<) 

1  F0RMATII3) 

3  FORMAT ( F7.1 ,ER.  1  ) 
10  FORMAT!1*'.      ROI'*')l 
22  FORMAT!'*1,  20X,1  (  INPUT  nATA  ) 

2S  FORMAT! ' 1 '. 10X , 'THIS  CASE  FOR 

?0  FORMAT! 'l'.lOX, 'THIS  f.ASF  FOR 

3*  FORMAT! 'I' ,IOX, 'THIS  CiSC  FOR 

50  FORMAT! • 1 ' ,inx. -TH!S  CSSF  FOR 

51  FORMAT ( >1 ', 10X, 'THI S  CASE  FOR 
40  FORMAT!  'I'.IOX,  'THIS  CASF  FOR 

52  FORMAT! >1' ,10X, 'THIS  CASF  FOR 
54  FORMAT! ' 1 ', 10X, 'THIS  CASE  FOR 

44  e«M«tl!*».  - .F12.5.5X..CFMS 

,Et2.5.5X.'STl 

.F12.5.5X, 'FE1 
.F12.5.5X,  '  SEFI. 
.F12.5.5X. >€ ICF 
,F12.5,5X,1OET 
.F12.5.5X, •OCHT 
.E12.5) 

•E12.5.5X. 'T4 
.E12.5.5X, 'AVE 
.F12.5.5X, 'GRO 
•F12.5.5X. 'APD1 
•  F12.5.5X,  'RR 
,C12.5,5X.'GR1 


200  FORMAT! 


102  FORMAT! 


104  FORMAT! 


106  FORMAT! 


/.OX, ' H 
/.ox, 'WT1 
/.OX, 'ST2 
/.OX  ,  'WEFI, 
/.OX, 'E  ICL 
/,0X, 'OGT 
/,0X, 'DPTR 
/.OX, 'V 
/.OX, 'T3 
/,QX, ' AVS 
/,0X, 'VS2 
/ ,ox, •ORE 
/,QX, ■RA 
/.OX,  '0E2 
/.OX, 'CFMS 
/,°X,'T3 
/.OX, 'AVE 
/.OX,  'VS2 
/,0X, "RR2 
/.OX, 'RA 
/.OX, '0E2 
/ ,QX, ■CFMF 
/,0X, 'TGRO 
/.OX, 'TOE 


AIR  VELOCITY  (V)=100  FPM  ' 

AIR  VELOCITY  <V)=200  FPM  ' 

AIR  VELOCITY  (V)=300  FPM  ■ 

AIR  VELOCITY  (V)=400  FPM  ' 

AIR  VEI0C1TY  IVI=500  FPM  ■ 

AIR  VELOCITY  (V)=fcOO  FPM  • 

AIR  VEIOCITY  IV)=700  rPM  ' 

AIR  VELOCITY  (V)=ROO  FPM  ■ 

' ,E12.5,5X,'CFME  = 
,E12.5,5X,'WT2   * 

'  .E12.5.5X, 'FE2  = 
,E1?.5,5X,'EICU  = 

'  .E12.5.5X, 'DOT  = 
,E12.5,5X,'DPKW  = 

'.E12.5.5X, 'OCCT  = 


*',E12.5,5X,'T5 
■' .E12.5.5X, 'VS1 
=  '  .F12.5.5X, 'ORG 
=  '  .E12.5.5X,' A°D2 
=',E! 2.5,5X,'0E1 
1  ,C1?.5,5X, 'GR2 


,Fil2.5.5X,  'CFMSl«'.E12.5t5X,'f  A 


E12.5.5X, 'T4 
.E12.5.5X, 'FA 
.F12.5.5X, 'GR 
■E12.5.5X 
,E12.5,5X 
■E12.5,5X,>6 
■E12.5) 
.  Fl 2.5,5 X, 'TGRG 


=  '  .E1P.5.5X,1 AVS   = 

=',E12.5.5X,'VS1   = 

=' ,E12.5,5X,'GR1   = 

APD1  =>,E12.5,5X,'APD2  = 

RB    =• ,E12.5.5X,'0E1   = 

G     =',E12.5,5X,'CFMS  = 


',E1?.5.5X,>T<",RE 


E12.5,5X.<TA0E  = '  ,E 12 .5 , 5X , ' ENRO 


,E12 

,E12, 
,E12 
,E12, 
,E12 
,E12 
,E12 


,E12. 
,E12. 
,E12. 
,E12. 
,E12. 
•  E12. 
,E12. 
,E12. 
,E12, 
,E12. 
,E12. 
,E12. 
,E12. 

,E12, 
,E12, 
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*  /,<)X,  'PAYRE= 

MM=MM+1 

IF(MM-KK)  5,5,7 
7  CONTINUE 

c    tSL.*.»*««.««^ 

c 
c 


IF(INnEX,EO,2)  GO  TO  31 

TMIN=WT1 

TMAX=WT3 

RO  TO  32 

31  TMIN=ST3 
TMAX=ST1 

32  CONTINUE 
NK=N-l 

DO    200    1=1. NK 

"if'^THIN   RE    0LnTn).4ND.TM!N.LE.ni.OTU+l)    J   CO   TO    1« 
IF     [TMAX!(;I;0I  DT(J-1).AND.TMAX.LE.0I.DT(J)  )     GO    TO    15 
GO    TO    200 
l*    ^,IN=(niDHIKKl  +  l»-OLnH(KKl))/(OLOT(KKltl)-nLOT(KKl))* 

I     <TMIN-0\0TIKK1))+0I_DH(KK1> 
15    55»X»rotnH(KK2*l)-Ol.nH(KK2))/IOLBT(KK2+l)-OI.DT(KK2))* 
1     (TMAX-ni_0T(KK2)  1+DI.0HIKK2I 
200    CONTINUE 

00    100    1=1. NK 

jflS?J|Xr:l?0?0HU)I2MDlHMIN...E.nU0H(I+X))    GO   TO   40 

iFIHMAxlGElniDHU-n.ANO.HMAX.LE.OLOHlJ))    GO    TO    30 

GO    TO    100 

A1    fF?HM?N,!l.B.niOH(I).»ND.HMIN.Re.nt.n^HU  +  l)J    GO    TO    40 
IFIHHAX.LE.OLOHIJ-ll.ANn.HMAX.GE.OLDHIJ))    GO    TO    30 

GO    TO    100 
40    K1=I 
30    <2=J-l 
100    CONTINUE 

NN=K2-K1 

NN1=NN+1 

NN2=NN+2 

RNFWH( 1 |=HMIN 

UNFWHfNN?) =HMAX 

RNEWTI 1 )=THIN 

UNFWT (MN? )=TMAX 

00    80    I_  =  2.NN1 

UNEWHIL  )=0I.0H(K1  +  L-1) 
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RNEWT(L)=OLDT(Kl+l-l> 
80    CONTINUE 

^]LTroRir^Siri*FTSfJ*'?THAX-RNEWTIL))*RNEWH<L) 
WT1=RNEWT(I.  +  1) 
WT10I.0=WT1 
T4HGR=T4HGR+HGR 
l.l.  =  LL+l 
GO    TO    17 
"5    rc,k'l    nn*rPM5*F«;    *(RNEUTINN2-J>-TMIN)*RNEWH<NN2-J) 
STl=nLDT(N'M4X-LLL) 
ST10i.O  =  5Tl 
TACGR=T4CGR+CGR 
I  I  I  =1  I.L+1 
97    CONTINUE 
RETURN 
END 
SENTRY 
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ABSTRACT 

Energy  conservation  is  now  an  integral  part  of  good  planning,  for 
efficient  use,  and  systems  are  being  recognized  which  are  economically 
and  technically  feasible.  Energy  conservation  may  therefore  become  compul- 
sory. 

Heat  lost  by  natural  ventilation  cannot  be  recovered,  and  the  only 
energy  saving  possibility  is  to  reduce  the  rate  to  the  unnecessary  wastage 
due  to  structural  design,  such  as  the  use  of  windows  with  controlled 
crackage. 

For  mechanically-ventilated  buildings  however,  using  inlet  and  exhaust 
systems,  a  number  of  techniques  exist.  It  is  possible  to  save  at  least 
half  of  the  heat  lost  in  the  exhaust  air.  These  techniques  generally  in- 
volve use  of  heat  recovery  devices  of  some  type. 

Heat  recovery  systems  can  be  grouped  basically  into  direct,  and  in- 
direct heat  transfer  apparatus.  In  the  direct  system,  the  warm  exhaust 
air  and  the  incoming  cold  air  pass  through  the  same  heat  exchanger  unit. 
The  transfer  between  the  two  air  streams  are  either  by  conduction  or  by 
regenerative  process  as  in  the  thermal  wheel. 

In  the  indirect  system,  a  secondary  medium  is  used  to  transfer  heat 
from  the  exhaust  and  cold  air  streams.  This  has  the  advantage  that  the 
two  air  streams  do  not  have  to  be  brought  together. 

Two  systems  which  are  typical  of  other  types,  run-around  coil  system 
and  a  heat-pipe  recovery  unit  were  selected  for  study.  These  were  selected 
for  analysis  because  of  their  highest  net  energy  recovered  with  respect  to 
a  reasonable  heat  recovery  unit  size  which  involves  optional  initial  cost. 


Effects  of  different  parameters,  glylcol  ethyene  and  air  flow  ve- 
locities for  the  run-around  coil  system,  and  the  effect  of  change  in  air 
flow  velocities  for  heat-pipe  recovery  unit  were  studied  for  highest  net 
energy  saved  and  pay  back  time. 

It  was  noticed  that  the  energy  recovered  increases  as  the  heat 
transfer  area  increases,  thus  increasing  initial  cost.  This  caused  an 
increase  in  the  operating  energy  cost  affecting  net  energy  saved.  The 
net  energy  saved  reaches  a  maximum  limit  with  respect  to  the  parameter 
affecting  the  change  in  the  heat  transfer  area  and  then  tends  towards  a 
decrease. 

Energy  saved  and  pay  back  time  depends  on  the  initial  costs  and  on 
heating  and  cooling  costs.  For  this  study  these  costs  were  selected  at 
current  levels  and  may  be  altered  in  order  to  assess  the  cost  involved 
for  a  localized  area. 

Keeping  other  parameters  constant  but  varying  only  one  at  the  same 
time  resulted  in  the  highest  net  energy  saving  for  a  specific  value  of 
that  parameter.  The  same  result  was  obtained  for  a  range  of  velocity  with 
respect  to  the  selected  values  and  the  local  design  temperature. 


